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PREFACE 


On  30  September  1992,  the  U.S.  Army  Ballistic  Research  Laboratory  was  deactivated  and 
subsequently  became  a  part  of  the  U.S.  Army  Research  Laboratory  (ARL)  on  1  October  1992. 
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1.  INTRODUCTION 


A  basic  mission  of  the  Advanced  Propulsion  Technology  Branch  (APTB),  Propulsion  and  Flight 
Division  (PFD),  U.S.  Anny  Research  Laboratory  (ARL),  is  the  evaluation  of  novel  propulsion  concepts 
and  propellant  formulations  which  often  involve  unusual  chemistries.  Given  the  unique  nature  of  many 
of  the  proposed  propulsion  concepts  and  the  lack  of  detailed  information,  especially  burning  rate, 
associated  with  many  of  the  candidate  propellants,  detailed  computer  simulation  is  often  impossible. 
However,  regardless  of  the  concept,  if  projectile  acceleration  is  achieved  through  a  gas  dynamic 
(thermodynamic)  process,  limiting  performance  estimates  can  be  obtained  through  the  use  of  a  constant 
breech  pressure  (CBP)  simulation. 

CBP  simulation  is  particularly  useful  in  quantifying  potential  performance  of  new  propellant 
fomiulaiions.  As  the  propellant  formulaior  explores  chemistries  beyond  the  traditional  single/doubleAtiple 
ba^  chemistry,  impetus  as  a  measure  of  potciuial  perfonnance  for  a  propellant  becomes  less  reliable. 
Only  a  ballistic  simulation  can  account  for  propellant  tliemiCKdiemical  properties  (i.e..  covolume,  raiio-of- 
specific  heats  (gamma),  etc.)  not  incorporated  in  the  impetus  calculation  but  which  caat  vary  over  a  much 
larger  range  than  associated  with  standard  propellant  formulations.  Iw  example,  the  value  of  gamma  for 
traditional  solid  prot>ellams  gcttcrally  varies  between  1.22  and  1.25.  yei  recently,  several  propcllaiu 
fommlaiiomi  rc?ve  been  proposed  wiUi  gammas  in  Ure  range  of  1 .05  to  1 . 1 . 

Ballisiically.  die  CBP  simulation  offers  several  advantages  in  cvaluatirrg  potential  gun  performance. 
First,  the  constsni  brccdi  pressure  gun  gives  Uie  maximum  possible  velodiy  whiclr  can  be  adheved 
without  lisc  usr  nr  'ectiniqucs  such  as  traveling  diargc  to  alter  dow-m  tube  pressures.  Tlius,  tire  calculation 
provides  an  absolute  m-^sure  of  maximum  vcltKity  perfonnance.  Next,  for  solid  propdlam  simulations, 
velocity  predictions  arc  dc{Kndeni  upon  Utc  grain  geometry  .selected  for  live  simulaiimi:  assuming  a  CBP 
eliminates  variaWc.  Finally,  trew  propulsion  co  cpts  often  involve  liquid,  gel  (liquid/Uiickcncr), 
emulsified  (liquiciyiiquid),  or  slurry  (.solid/liquid)  propellants  which  can  te  evaluated  only  by  a  CBP 
simulation  due  to  tire  lack  of  fixed  propellant  geometry.  In  addition,  it  has  been  .shown  (Irisli  1985)  Uiat 
well  designed  solid  piojreUani  gun  systems  achieve  velocities  between  90%  and  95%  of  the  CBP  velocity. 
Tabic  I  illustrates  ihi.s  fact  for  a  variety  of  fielded  gun  systems  (Morrison  1990).  Tire  ratio  of 
experimental  vclotrity  to  CBP  velocity  is  termed  live  ’'ballistic  ratio”  (BR). 
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Table  1.  Experimental  Velocity  and  Constant  Pressure  Velocity  for  Various  Cannons 


Gun 

Caliber 

(mm) 

Propellant 

Type 

Chamber 
Vol.  (1) 

Pmax 

(MPa) 

Vel  Exp 
(m/s) 

Vel  CP 
(m/s) 

Bofors  L70 

40 

NC1066 

0.55 

319 

1,005 

1,079 

0.93 

IMl 

60 

M30 

2.33 

460 

1,620 

1,781 

0.91 

M68 

105 

M30 

6.47 

414 

1,486 

1,620 

0.92 

M256 

120 

JA2 

9.75 

510 

1,650 

1,740 

0.95 

XM25 

120 

JA2 

9.75 

510 

1,739 

1,828 

0.95 

XM25* 

120 

JA2 

10.00 

683 

2,423 

2,490 

0.97 

Navy  5754 

127 

NACO 

13.(W4 

372 

808 

860 

0.94 

M198 

155 

M30A1 

I'lji 

18.85 

313 

826 

884 

0.93 

Hig))  v«locUy  test  with  kg  humch  mus 


The  objective  of  litis  report  is  to  documem  tlte  theory  and  operation  of  Ute  constant  breech  pressure 
computer  code  CONPRESS.  In  addition  to  ilie  standard  CBP  calculation,  CONPRESS  incorporates  a 
number  of  optimization  and  search  algoriUuns  together  with  extensive  parametric  variation  capabilities. 
Documentation  and  use  of  tltesc  t^ons  is  also  provided. 

2.  THEORY  AND  DERIVAHON  OF  EQUATIONS 

(NOTE:  Equivalent  derivations  have  been  pcrfomied  by  other  rescarcliers  over  the  paa  50  years  and 
can  be  found  in  various  pubheattorvs.) 

a  0.  K  Conaanis 

A  Projectile  base  area 

Cp  Heat  capacity  at  consiaiu  pressure 

Heat  capacity  at  cotrstani  volume 
e  Imcntal  gas  energy 

F  Force  acting  on  projectile  base 

7  Ratio  of  specific  heats 
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M 

m, 

P 


P 


P. 


m 


be 


R 

S 

T 


U, 


m 


'm 

n/ 

V;./ 

Xj 

Xb 

*tn 


Propellant  impetus 
Kinetic  energy  of  projectile 
Kinetic  energy  of  gas 
Molecular  weight  of  gas 
Chaige  mass 
Projectile  mass 
Propellant  covolume 
Chamber/breech  pressure 
Constant  breech  pressure 

Space  mean  pressure 

Constant  space  mean  pressure  during  propellant  burning 

Space  mean  pressure  at  projectile  exit 
Base  pressure 

Constant  base  pressure  during  propellant  burning 

Universal  gas  constant 

Entropy 

Gas  temperature 

Gas  temperature  at  burnout 

Flame  temperature  of  propellant 

Projectile  muzzle  velocity 

Total  volume  as  a  function  of  projectile  position 

Gun  volume  from  0  to  x^,,  chamber  and  tube  volume  to  propellant  burnout  position 

Gun  volume  from  0  to  Xj,  chamber  volume 

Gun  volume  from  0  to  x„,  chamber  and  tube  volume 

Free  gun  volume  at  burnout  adjusting  for  propellant  covolumc 

Free  chamber  volume  adjusting  for  propellant  covolumc 

Free  gun  volume  at  projectile  muzzle  exit  adjusting  for  propellant  covolume 

Initial  projectile  position 

Position  of  projectile  at  propellant  burnout 

Projectile  position  at  muzzle  exit 


For  the  derivation,  the  following  assumptions  arc  assumed: 
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(1)  Tile  Lagrange  gradient  adequately  describes  the  pressure  and  gas  velocity  profiles  of  the  ballistic 
cycle.  Equations  1-3  summarize  the  relations  due  to  Lagrange  v.hich  will  be  used  in  the 
derivation. 


m 

1  +  _i 

2ot, 


Ph. 


1  +  — L 

V 


K  ^  K 
«  Tni; 


(1) 


(2) 


(3) 


(2)  The  propellant  bums  in  an  ideal  manner  (i.e..  the  p.  opellant  is  instantaneously  converted  to  a 
gaseous  state). 

(3)  The  burn  rate  of  the  propellant  can  be  controlled  so  as  to  provide  a  constant  chamber  pressure  till 
propellant  burnout, 

(4)  The  gas  is  polyuoplc, 

(5)  The  Nobcl-Abcl  equation  of  state.  Equation  4.  is  valid. 

(4) 

(6)  No  energy  losses  occur  during  die  ballistic  cycle. 

(7)  After  burnout,  adiabatic  expansion  of  the  gas  takes  place. 

(8)  Hie  projectile  base  area  equals  the  cross>secUonal  area  of  the  tube. 
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A  diagram  of  the  gun  with  accompanying  coordinate  system  is  given  in  Figure  1.  The  location,  x^, 
conesponds  to  tube  entrance,  to  the  location  of  the  projectile  at  propellant  burnout,  and  to  muzzle 
exit 


End  of  Chamber 


- - ^ 

Chamber 

1  t 

0  X, 


Tube  Burnout  Muzzle 

Entrance  Poeltion  Exit 

Figure  1.  Coordinate  system  for  (he  gun. 

Under  the  assumption  of  no  energy  losses,  ite  kinetic  energy  of  (he  i»ojcctile  as  it  travels  from  its 
initial  position,  x^,  to  muzzle  exit,  is  given  by  Equation  S. 


a:^ (5) 

where  F(x)  is  the  force  on  Uic  projectile  base  at  position  x,  x,  S  x  S  x^.  Note  that  Equation  5  is  not 
dependent  on  tltc  CBP  assumption  but  is  simply  the  deOnition  of  wotk.  If  A  is  the  projectile  base  area 
and  Pfjtx)  base  pressure  at  position  x,  Xj  x  £  x„,  then 


F(x)->U>^(x)  . 


(6) 
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Also,  fiom  the  assumption  that  the  projectile  base  area  is  equal  to  the  cross-sectional  area  of  the  tube. 


V(jc)  -  A(x  -  x^)  +  V^=>dV  -  Adx  , 


(7) 


where  V(x)  is  the  total  volume  of  chamber  plus  gun  tube  to  location  x,  xi,ix^  x^.  Substitution  of 
Equations  6  and  7  iiito  Equation  S  yields 

AP^(x)dx^  j^^P^(V)dV  ,  (8) 


where  =  V(x^)  and  =  V(xJ. 

If  burnout  occurs  at  or  after  muzzle  exit,  then  the  integral  in  Equation  8  can  easily  be  evaluated  if  a 
CBP  is  assumed  since  by  the  Lagrange  a<:^mption.  Equation  1.  the  base  pressure  also  remains  <»nstmiL 
This  gives 

->l)v  <« 

or,  using  Equation  1  to  truisfomt 
be  given  by  Equation  10, 


However,  if  burnout  occurs  before  muzzle  exit,  Ute  ba<%  pressure  will  have  a  pfofilc  as  sltown  in 
Figure  2.  and  tltt  energy  impaiicd  to  the  projectile  during  >w(pansion  piuise.  'rom  to  x^,  must  also 
be  calculated.  From  the  assumptions,  it  is  assumed  that  the  crnLision  will  be  adiabatic.  Thus,  the  work 
performed  on  the  projectile  by  the  expanding  gas  is  a  fbmoion  of  sp^re  mean  pressure  instead  of  base  or 
chvnbcr  pressure.  Thus.  Equation  8,  which  giv&<‘  the  u>tai  kiiselic  energy  of  the  projectile  in  terms  of  base 
pres^re,  must  be  rewritten  in  terms  of  ^»cc  means  pressure.  To  puform  this  transformation,  consider 
multiplying  both  sides  of  Equation  8  by  the  factor  (1 


from  base  to  breech  pressure,  the  kinetic  energy  of  the  projectile  would 


^  ( .*1(1  -  ) 


m, 

2. .. 


(10) 
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iroieciiie  base  nfessune  vs.  piotecu'te  travel  if  propellant  burnout  occurs  prior  to  muzzle  exit 


(11) 


/ 

m. 

1  +  -i- 

- 

3m 

P 

V 

A 


Using  Equation  2  yields. 


f 

1  + 


jl”  P{V)dV. 


(12) 


It  is  of  interest  to  note  that  the  integral  on  the  right-hand  side  of  Equation  12  leprrsents  the  total 
usable  work  which  can  be  extracted  from  the  propellant,  while  the  left-hand  side  is  the  sum  of  the  kinetic 
energy  of  the  projectile  and  the  kinetic  energy  of  the  gas  using  the  Lagrange  assumption  given  by 
Equation  3  to  relate  the  kinetic  energy  of  the  gas  and  projectile. 

From  tliis  point  forward,  a  constant  breech  pressure  during  propellant  burning  is  assumed. 
Equation  12  provides  a  means  for  determining  Kp  by  evaluating  the  integral  on  the  right-hand  side  of  the 
equation.  While  the  propellant  bums,  by  the  Lagrange  assumption,  the  chamber,  space  mean,  and  base 
pressures  remain  constant.  Letting  be  the  constant  space  mean  during  the  time  the  propcUam  bums, 
Equation  12  can  be  rewritten  in  the  following  manner. 


f 

1  > 


m 


c 


J^-F(V)dv. 


Evaluating  the  first  integral  gives 


/ 

1 


(13) 


(14) 


To  evaluate  the  second  integral,  P(  V )  must  be  replaced  with  an  expression  involving  V  as  the  variable. 
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For  a  polytropic  gas  (Whitman  1974)  during  the  expansion  phase, 

s 

T  (15) 

PiV-i\m^)^  »  ae  '  , 

where  y  is  the  specific  heat  ratio,  ri  the  propellant  covolume,  S  the  entropy,  the  value  of  specific  heat 
at  constant  volume,  and  a  a  constant.  Also,  for  an  adiabatic  gas,  entropy  is  constant  Thus,  assuming 
remains  constant  Equation  IS  becomes 

F(V-ilw^)^-/:,  (16) 


where  K  is  a  constant. 


Solving  for  P  yields 


P  - 


K 


Substituting  Equation  17  into  Equation  14, 

(V-nm/ 


(  \ 

\  4 

3m,  ^ 


(17) 


(18) 


Integrating, 


I  , 

Jm, 


(y-l)(V'-n/n,)Y’i  ’ 


(19) 


or 


1  ♦ 


JUi 

3m 


P  J 


(y-l)(l'^-nmc)^“*  (Y-  l)(V'^-Tlm^)l^"' 


.  (20) 


To  simplify  Equation  20,  the  last  two  terms  can  be  rewriuen  to  give 
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(y-i)iv^-nm^)'^  (y-i)(v^-nw,)Y 


.  (21) 


which  upon  substitution  of  Equation  17  gives  Equation  22, 


1  *  3^  k'-Pc(^6-''c)--^-^ - -y-r  "•• 


By  assumption,  the  constant  chamber  pressure  duririg  propellant  burning  is  known  and  thus  is  also 
known  by  the  Lagrange  assumptions.  The  geometry  of  the  gun  fixes  and  V^.  Finally,  covolume,  fi, 
gamma,  y,  and  propellant  mass,  m^,  are  known  properties  of  the  propellant,  and,  hence,  the  only  unknown 
quantities  in  Equation  22  are  Kp,  and  P^.  To  eliminate  P^,  apply  Equation  16  at  the  beginning  and 
end  of  the  expansion  phase.  This  gives 

»»c)^  •  ^  •h  (yi>  -  n«c  )^  . 


i  .i f y 

“  ‘  v--nm;  • 


Substitution  of  Equation  24  into  Equation  22  yields 


3'np]  '  Y“l  Y-1  "  ' 
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For  convenience,  let  the  frae  volumes  be  given  by 


Vcf  V'c  -nw,, 


Vbf  Vb 


(29) 

Finally,  solving  Equations  1  and  2  for  space  mean  pressure  in  terms  of  chamber  pressure,  and  applying 
tire  result  to  the  kiwwn  constant  chamber  pres^ire,  produces  Equation  30, 


1 

3m„ 


>.J!i 

2m, 
\  f 


Substituting  Equations  26-30  into  Equation  25, 


I  *  ^  - 

3m  J  ^ 


V 

“KT/ 


which  upon  simplifying,  gives  an  expressing  for  K.  in  Equation  32. 
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(32) 


f 

Pc 

■<  V..-V.-  VI 

fvl 

Y 

1  +  —L 

1 

1 

1 

> 

Intj, 


Before  determining  an  expression  for  V  the  consistency  of  the  equations  can  be  checked  by 
assuming  burnout  at  muzzle  exit  (i.e.,  in  Equation  32).  This  should  reduce  Equation  32  to 

Equation  10.  Letting  =  V^^in  Equation  32  yields 


V  V 

Vmf-  ^cf-  Yn 


V 


Vj 


(33) 


Sim{difying,  Equation  10  is  obtained: 


(34) 


To  obtain  an  expression  for  Uic  Nobcl'Abcl  equation  of  state.  Equation  4,  together  with  energy 
conservation  during  the  time  which  the  propellant  is  iKiming.  will  be  utilized.  At  bumouu  the  Nobcl'Abcl 
equation  of  stats  becomes 


m. 


fin 

M 


(35) 


Note  that  the  unknown  gas  temperatute,  has  been  imroduced  and  needs  to  be  eliminated.  Now. 
while  the  propellam  bums,  since  no  energy  loss  is  assumed,  the  change  in  the  imcmal  gas  energy  must 
translate  into  the  work  being  perfomicd  on  the  projectile  and  propelling  gas.  The  dtangc  in  internal 
energy  of  a  gas  per  unit  mass  is  given  by  Equation  36  (Whitman  1974). 
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(36) 


e  =  c^dT . 

Thus,  the  total  change  in  internal  energy  during  propellant  burning,  assuming  a  constant  c^,,  is 

mf.t  =‘ m^c^dT  «  m^c^,(7y- r^)  .  (37) 

Now  the  work  done  on  the  gas  and  projectile  during  the  burning  of  the  propellant  is 


^c(  W- 


(3«) 


Equating  the  expressions  for  work  in  Equation  37  and  38  gives 


^(n/- V'”cfv(7>-n). 


(39) 


Equations  35  and  39  provide  a  system  of  two  equations  in  two  unknowns.  W^^tmd  which  can  be  solved 
for  as  given  in  Equation  40. 


m. 


M 


^  M  c  cf 


(40) 


By  dcTmiiion.  propellant  impetus  or  force,  /,  Is  given  by 


/  - 


(41) 


and 
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(42) 


A 

Af 


Substituting  Equations  41  and  42  into  Equation  40  and  simplifying. 


'bf 


mgCj  ^  (Cp-cJV^f 


or 


(43) 


(44) 


Making  use  of  Equation  30  to  eliminate  the  final  equation  for  is  obtained: 


m^l 


Rt. 

1  ♦ 

2  m 


1  ♦ 


JUl 

3  m, 


L.  .  IJ.  V 

^  Y 


(45) 


Equation  32  together  with  Equation  45  provides  an  expression  for  in  terms  of  P^.  m^,  m^,  y.  /.  ^ef 
V  However,  and  can  both  be  expressed  in  terms  of  and  arc  given  in  Equation  46  and  47: 


(46) 

(47) 


Thus,  Equations  32  and  45-47  provide  a  means  to  compute  in  terms  of  the  known  quantities  P^. 
m^,  ntp,  Y.  /.  A,  V^.  q,  and  -  Xj).  the  projectile  travel.  The  muzzle  velocity  can  then  be  computed 
using 
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It  is  also  of  interest  to  know  the  location  of  the  projectile  at  propellant  burnout.  From  Equation  28, 


Substituting  in  Equation  45  and  solving  for  j. 


Equation  50  can  be  used  to  detemtine  the  projectile  position  at  propellant  bum-out  relative  to  the  inid^ 
position  of  the  projectile. 

3.  PROGRAM  OPTIONS,  REQUIRED  INFORMATION.  AND  INPUl'  DATA  PIECK  STRUCTURE 

Tlic  equations  derived  iit  the  previous  section  together  with  optimization  algonthms  have  been  cod^ 
in  FORl’RAN  to  produce  the  computer  code  CONPRESS  (Appendix  A).  Seven  different  computational 
options,  described  in  TaUc  2.  are  available  in  CONPRESS  Table  3  summarizes  the  required  input 
infonnaiion  for  each  option.  For  option  6  tltc  required  input  is  determined  by  the  option  (I -5)  sdeoed 
for  the  pantm^iic  study. 
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Table  2.  Computational  Options  for  COMPRESS 


Computational 

-  ....1  ,  iW-  IJ-  .  ^B-  ■  _L  1  LM 

Description 

1 

Compute  muzzle  velocity  given  chamber  volume  and 
charge  mass. 

2 

Compute  required  charge  mass  given  chamber  volume, 
desired  velocity,  and  maximum  loading  density. 

3 

Compute  required  charge  mass  and  chamber  volume 
given  desired  velocity  and  loading  density. 

4 

Compute  the  charge  mass  to  produce  optimal  muzzle 
velocity  given  the  chamber  volume  and  maximum 
loading  density. 

5 

Compute  the  charge  mass  and  chamber  volume  to 
produce  optimal  muzzle  velocity  given  loading  density. 

6 

Performs  a  parametric  study  over  a  specified  range  of 
one  or  more  of  the  input  variables,  the  study  can  be 
perfomied  for  any  of  the  options  1-5. 

7 

Option  to  end  program. 

Tabic  3.  Required  Input  Data  for  Computational  Options  1-5 


The  structure  of  the  input  deck  for  CONPhiiSS  for  single  and  muldple  runs  is  given  in  Table  4. 
[Note:  The  commas  utilized  in  Table  4  to  sepamte  entries  on  the  same  line  are  not  required  in  an  actual 
input  deck,  a  space  should  separate  fh.  entries.]  It  should  also  be  noted  that  all  data  items  must  be 
assigned  a  value  even  if  that  parameter  is  not  used  in  the  computation.  For  example,  option  1  does  not 
require  input  of  the  loading  density  since  both  ht;  charge  mass  and  chamber  volume  are  input  values. 
However,  a  value  must  be  entered  in  the  input  deck  for  the  loading  density.  This  value  will  be  ignored 
in  the  computation.  In  the  following  sections  additional  infomiation  about  specifics  of  the  input  deck  for 
each  option  will  be  discussed. 

As  indicated  in  Table  4,  if  multiple  runs  are  to  be  perfonned,  an  output  file  name  is  not  repeated  for 
each  run.  All  output  is  written  to  the  same  file.  Note  also  that  to  terminate  the  program,  a  7  should  be 
entered  after  the  last  data  deck,  no  title  is  needed. 


Table  4.  Structure  of  Input  Deck  for  CONPRESS 


Single  Run 

Line  1:  Output  file  name 

Line  2:  Title 

Line  3;  Run  Option  (st  Table  2) 

Line  4:  Maximum  Breech  Pressure  (Constant  Breech  Pressure  Value) 

Line  5:  Charge  ‘  lass.  Projectile  Mass 

Line  6:  Gamma,  Covolume,  Impetus,  Propellant  Density 

Line  7:  Chamber  Volume,  Bore  Diameter,  Projectile  Travel 

Line  8:  Desired  Velocity,  Loading  Density 

Line  9:  7  (to  indi''ate  end  of  run,  program  will  correctly  terminate  if  this  line  is  not  included. 

Multiple  Runs 

Line  1:  Output  file  name 

Line  2:  Title  Run  1 

Line  3:  Run  Option  (see  Table  2) 

Line  4:  Maximum  Breech  i  ressure  (Constant  Breech  Pressure  Value) 

Line  5:  Charge  Mass,  Projectile  Mass 

Line  6:  Gamma,  Covolume,  Impetus,  Propellant  Ovjnsity 

Line  7:  Chamber  Volume,  Bore  Diameter,  Projectile  Travel 

Line  8:  Desired  Velocity,  Loading  Density 

Line  9:  Title  Run  2 

Line  10:  Run  Option  (see  Table  2) 

Line  11:  Maximum  Breech  Pressure  (Constant  Breech  Pressure  Value) 

Line  12:  Charge  Mass,  Projectile  Mass 

Line  13:  Gamma,  Covolume,  Impetus,  Propellant  I>cn.':iiy 

Line  14:  Ciiamber  Volume.  Bore  Diameter,  Projectile  Travel 

Line  15:  Desired  Velocity,  Loading  Density 

Lines  9-15  are  repeated  for  each  subsequent  mn 

Line  N:  7  (to  indicate  ci»d  of  run,  program  will  correctly  terminate  if  this  line  is  not 

included). 
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TESl'  CASE  1:  Computational  Option  1 


Option  1  is  a  direct  application  of  the  equations  to  predict  projectile  velocity,  kinetic  energy  and 
projectile  position  at  propellant  burnout.  Both  entries  on  line  8  (desired  velocity  and  loading  density)  of 
the  input  file  are  not  utilized  in  the  computation.  However,  as  mentioned  in  the  previous  section,  values 
for  each  must  be  provided.  Table  5  provides  a  sample  input  deck  for  option  1,  the  corresponding  output 
for  the  COMPRESS  run  is  given  in  Table  6. 


Table  5.  COMPRESS  Input  Deck  for  Option  1 


cpv21.otl 

'H'is  tests  optiori  1  for  latest  version  of  COMPRESS. 
1 

575. 

9.  9. 

1.225  .990  1140.  1.58 
10.5  120.  475. 

1200.  .95 
7 

tq  ■  iiiisnxi  M  m  ;.iw.  I  ■;  'i  i  n  ■  gastw;  as;  j  utp.  -  incases 


To  verify  the  coding  the  interior  ballistic  code,  IBHVG2  (Anderson  ijid  Fickic  1987)  was  run  with 
the  same  input  parameters  (Table  5)  using  the  cottsiant  breech  pressure  mode  with  a  Lagrange  gradient. 
Pertinent  re  sults  of  the  IBH  VG2  nrn  ate  given  in  Table  7.  As  crat  be  seen  from  a  comparison  of  Tables 
6  and  7  results  are  virtually  ide.ntical. 

TEST  CASE  2:  Computational  Option  2 

In  option  2,  the  charge  mass  .  ^tquired  to  obtained  tlic  desired  velocity  entered  on  line  8  of  die  input 
deck  is  the  quantity  sought.  Thus,  all  input  values  except  for  tire  cliarge  mass  enteral  on  line  5  arc 
utilized.  A  simple  bisection  mcUrod  for  root  finding  is  cmployal  to  obtain  tire  desired  charge  mass.  The 
function  used  is  tire  dis.crcncc  between  tite  pfeiicied  velocity  and  the  desired  velocity  as  a  frmetion  of 
charge  mass.  The  initial  interval  for  titc  diargc  mass  over  which  tire  search  will  be  perfonned  is  that 
obtairred  by  assuming  Urat  the  desired  velocity  will  occur  for  a  charge  mass  between  10%  and  100%  of 


Table  6.  COMPRESS  Output  for  Input  Deck  of  Table  5 


INPUT  PARAMETERS 

Input  File:  cpv21.opl 

Run  Option;  1 

Output  File;  cpv21.otl 

Gun  Parameters 

Prooellant  Parameters 

Cham.  Vol.  (L):  10.5000 

Impetus  (J/g):  1,140.00 

Tube  Diam.  (mm):  120.0000 

Gamma  (-):  1.22500 

Proj.  Tiav.  (cm):  475.0000 

Covolume  (cm^/g):  .99600 

Proj.  Mass  ^g):  9.00000 

Density  (g/cm^):  1.58000 

Max  Press.  (MPa):  575.0000 

Input  Prop,  (kg);  9.00000 

Tatget  Vel.  (m/s):  UOO.OOOO 

Target  LD  (g/cm^);  .95000 

RESULTS 

Ballistic  Values 

Energy  Information 

Velocity  (m/s);  1.667.5470 

Prop,  Energy  (MJ):  45.59999 

Burnout  (cm);  133.8050 

Proj.  KE  (MJ);  12.51320  (27.4%) 

1  Cham.  Vol.  (L):  10.50000 

Gas  KE  (MJ):  4.17107  (9.1%) 

Prop.  Mass  (kg):  9.00000 

Prop,  LD  (g/cm^);  0.85714 
Expansion  Ratio  (•):  6.11631 
Muzzle  Exit  Press  (MPa):  88.31 

Gas  Internal  (MJ):  28.91572  (63.4%) 

Mass  Chaige/Mass  Proj  (cAn  ratio): 
Piezomctric  Efficiency:  0.405 

1.00000 

Table  7.  RcsulUi  From  IBHVG2  Run  With  Identical  Inputs  for  Constant  Pressure  Calculations 


Ballistic  Values 


Velocity  (m/s);  1.667.67 
Burnout  (cm):  133.8 


Prop.  LD  (g/cm^):  0.857143 

Expansion  Ratio  (-)^  6.1 16 

Muzzle  Exit  Press  (MPa):  88.2791 

Mass  '  iiargc/Mass  Proj  (c/m  ratio):  1.00000 

Piczomeiric  Efficiency:  0.405 


Energy  Information 

Prop.  Energy  (MJ):  45.6 

Proj.  KE  (MJ):  1 2.5 1 5 1 20  (27.45%) 

Gas  KE  (MJ);  4.171707  (9.15%) 

Gas  Internal  (MJ):  28.907010(63.9%) 
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the  maximum  charge  mass  which  can  be  placed  in  the  chamber.  This  maximum  charge  mass  is  obtained 
by  multiplying  chamber  volume  (first  entry  line  7)  by  loading  density  (second  entry  line  8). 
Unfortunately,  this  choice  for  the  initial  charge  mass  interval  may  result  in  a  situation  where  the  charge 
masses  associated  with  endpoints  do  not  produce  a  sign  change  in  the  function  as  required  by  the  bisection 
method.  In  this  case,  the  user  is  warned  of  the  difficulty  and  is  prompted  for  a  new  interval  on  which  the 
search  is  to  be  perforaied.  This  condition  will  most  often  occur  when  the  desired  velocity  is  close  to  the 
optimal  velocity.  One  method  to  alleviate  this  problem  is  to  first  run  Option  4  to  determine  the  charge 
mass  for  optimal  velocity,  and  then  use  this  value  as  the  right-hand  endpoint  of  search  interval  when 
prompted  by  the  program. 

To  validate  the  search  algorithm,  1BHVG2  was  run  in  the  constant  pressure  mode  with  a  charge  mass 
of  7  kg.  All  other  input  remained  the  same  as  used  in  the  test  case  for  Option  1.  The  calculated  IBHVG2 
velocity  was  1,549.07  m/s  which  was  then  entered  as  the  desired  velocity  for  the  COMPRESS  calculation. 
The  COMPRESS  input  deck  is  given  in  Table  8  with  the  results  in  Table  9.  [Mote:  A  ran  time  error 
occurring  during  the  execution  of  the  program  indicates  that  the  desired  velocity  cannot  be  obtained  within 
the  design  parameters  (projectile  mass,  bore  diameter,  projectile  travel,  chamber  volume,  propellant 
properties,  and  maximum  pressure)  specified.] 

Table  8.  COMPRESS  Input  Deck  for  Option  2 
cpv2l.ol2 

I1iis  will  validate  option  2  of  latest  version  of  COMPRESS. 

2 

575. 

I.  9. 

1. 225  ,996  1140.  1.58 
10.5  120.  475. 

1549.07  .95 
7 


As  can  be  seen  in  Tabic  9,  COMPRESS  computed  a  charge  mass  of  7.00109  kg  to  obtain  the  desired 
velocity  of  1,549.07  m/s.  This  represents  a  0.015%  difference  with  the  IBHVG2  value  of  7  kg.  Thus, 
Option  2  appears  to  correctly  deicrmirtc  the  desired  charge  mass. 
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Table  9.  COMPRESS  Output  for  Input  Deck  of  Table  8 


INPUT  PARAMETERS 

Input  File:  cpv21.op2 

Run  Option:  2 

Output  File:  q)v21.ot2 

Gun  Parameters 

Prooellant  Parameters 

Cham.  Vol.  (L):  10.5000 

Impetus  (J/g):  1,140.00 

Tube  Diam.  (mm);  120.0000 

Gamma  (-):  1.22500 

ProJ.  Trav.  (cm):  475.0000 

Covolume  (cm^g);  0.99600 

Proj.  Mass  (kg):  9.00000 

Density  (g/cm^):  1.58(XX) 

Max  Press.  (MPa);  575.0000 

Input  ^p.  (kg):  1.00000 

Target  VeL  (m/s):  1349.0700 

Target  LD  (g/cm^);  0.95000 

RESULTS 

Ballistic  Values 

Encrev  Infonnation  I 

Velocity  (m/s):  1,549.0700 

Prop.  Energy  (MJ):  35.47216  | 

Burnout  (cm):  85.0453 

Proj.  KE  (MJ):  10.79828  (30.4%) 

Cham.  Vol.  (L);  10.50000 

GasKE(MJ):  2.79999(7.9%) 

Prop.  Mass  (kg):  7.00109 

Prop.  LD  (g/cm^):  0.66677 
Expansion  Ratio  (-).  6.11631 
Muzzle  Exit  Press  (MPa):  68.27 

Gas  Internal  (MJ):  21.87390(61.7%) 

Mass  Chaigc/Mass  Proj  (cAn  ratio); 
Piczomctric  Efficiency:  0.350 

0.77790  j 

TEST  CASE  3:  Computational  Option  3 

This  option  is  similar  to  Option  2  except  that  the  chamber  voliunc  is  also  a  variable.  The  calculation 
will  detennine  both  the  charge  mass  and  minimal  chamber  volume  <i.c.,  smallest  chamber  volume  which 
will  hold  the  cliarge  mass,  charge  mass  divided  by  loading  density)  required  to  obtain  the  desired  velocity. 
As  in  Option  2,  a  bisection  method  is  utilized  to  obtain  the  desired  charge  mass.  The  input  charge  mass, 
line  5,  is  taken  as  the  left-hand  endpoint  for  tlu;  initisd  search  interval.  The  right-hand  endpoint  is  obtained 
by  increasing  the  input  charge  mass  by  successive  multiplications  of  a  factor  of  1.2  until  a  suitable  value 
for  the  bisection  method  is  detennined.  The  chamber  volume  entered  on  line  7  is  not  utilized  in  the 
compulation.  [Note:  As  with  Option  2.  a  run  time  cnoi  occurring  during  the  execution  of  the  program 
indicates  that  the  desired  velocity  cannot  be  obtained  within  the  design  parameters  (projectile  mass,  bore 
diameter,  projcciilc  travel,  propcllam  properties,  and  maximum  pressure)  ^xcifled.] 
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Ir^ut  and  output  for  Option  3  are  given  in  Tables  10  and  1 1.  From  Table  1 1  the  desired  velocity  of 
1,500  m/s  is  obtained  using  5.39472  kg  of  propellant  with  a  chamber  volume  of  5.67865  1.  Running 
COMPRESS  with  Option  1  with  these  values  for  charge  mass  and  chamber  volume  yields  a  velocity  of 
1,500  m/s,  thus  validating  the  algorithm  for  Option  3. 

Table  10.  COMPRESS  Input  Deck  for  Option  3 
cpv21.ot3 

'liiis  is  the  test  of  option  3  of  latest  version  of  COMHiESS. 

3  I 

575.  I 

1.  9. 

1.225  .996  1140.  1.58 
10.5  120.  475. 

1500  .95 
7 


TEST  CASE  4:  Computational  Option  4 

Option  4  involves  the  dctcimination  of  charge  mass  which  will  result  in  optimal  muzzle  velocity  given 
all  other  input  parameters  fixed.  This  is  a  standard  mathematical  optimization  problem  (i.e.,  optimize  a 
function,  velocity,  of  a  single  independent  variable,  charge  mass).  Traditionally  two  approaches  are 
available  to  perform  an  optimization  of  a  function  of  a  single  variable.  One  involves  the  computation  of 
Utc  derivative  while  the  second  does  not.  For  COMPRESS,  the  former  approach  has  been  rejected  due 
to  the  complexity  of  taking  the  derivative  of  the  energy  equation  (Equation  32).  Instead,  an  inverse 
parabolic  interpolation  method  is  employed  in  searching  for  the  extrema.  Specifically,  the  method  utilized 
in  tile  program  is  a  method  due  to  Brent  (1973). 

In  the  input  deck,  both  the  charge  mass,  line  5,  and  the  desired  velocity,  line  8,  arc  ignored.  The 
interval  for  the  charge  mass  over  which  the  search  will  be  conducted  runs  from  0  kg  to  the  maximum 
mass  of  propellant  which  can  be  loaded  into  the  chamber,  chamber  volume  multiplied  by  loading  duisity. 

A  sample  input  deck  and  resulting  COMPRESS  output  arc  shown  in  Tables  12  and  13. 
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Table  11.  COMPRESS  Output  for  Input  Deck  of  Table  10 


1  IMPUT  PARAMETERS 

Ii^utFile:  cpv21.op3 

Run  Option:  3 

Output  File:  cpv21.ot3 

Gun  Parameters 

Prooellant  Parameters 

Cham.  Vol.  (L):  10.5000 

Impetus  (J/g):  1,140.(X) 

Tube  Diam.  (mm):  120.0000 

Gamma  (-):  1.22500 

Proj.  Trav.  (cm):  475.0000 

Covolume  (cm^/g):  0.99600 

Proj.  Mass  ^g):  9.0(XXX} 

Density  (g/cm^):  1.58000 

Max  Press.  (MPa):  575.0000 

Input  Prop,  (kg):  1.000(X) 

Target  Vel.  (m/s):  UOO.OOOO 

Target  LD  (g/cm^):  0.95000 

RESULTS 

Ballistic  Values 

Enercv  Information 

Velocity  (m/s):  1.500.0000 

Prop.  Energy  (MJ):  27.33324 

Burnout  (cm):  81.4230 

Proj.  KE  (MJ):  10.12500  (37.0%) 

Cham.  Vol.  a):  5.67865 

GasKE(MJ):  2.02302  (7.4%) 

Prop.  Mass  (kg):  5.39472 

Prop.  LD  (g/cm^):  0.95000 

Gas  Internal  (MJ):  15.18522(55.6%) 

Expansion  Ratio  (•):  10.46021 
Muzzle  Exit  Press  (MPa):  52.71 
Mass  Charge/Mass  Proj  (cAn  ratio): 
Piezometric  Efficiency:  0.328 

0.59941  1 

Table  12.  COMPRESS  Input  Deck  for  Option  4 
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Table  13.  COMPRESS  Output  for  Iiqjut  Deck  of  Table  12 


I 


For  a  chamber  volume  of  10.5  and  a  loading  density  of  0.95  g/cm^,  successive  nms  of  COMPRESS 
utilizing  Option  1  indicates  that  the  maximum  velocity  will  occur  when  the  chamber  is  totally  filled  with 
propellant.  Tliis  corresponds  to  9.975  kg  of  propcUam  (10.5  *  .95)  which  is  obtained  by  COMPRESS 
using  Option  4  (sec  Table  13).  Thus,  the  algorithm  for  Option  4  appears  to  be  functioning  correctly. 

TEST  CASE  5:  Computational  Option  5 


Input  File:  cpv21.op4 
Run  Option:  4 

Gun  Parameters 

Cham.  Vol.  (L):  10.5000 
Tube  Diam.  (mm):  120.0000 
Proj.  Trav.  (cm):  475.0000 
Proj.  Mass  (kg):  9.00000 
Max  Press.  (MPa):  575.0000 
Target  Vcl.  (m/s):  1.500.0000 


Balli.stic  Values 


IMPUT  PARAMETERS 

Output  Rle:  cpv21.ot4 

Propellant  Parameters 

Impetus  (J/g):  1.140.00 
Gamma  (-):  1.22500 
Covolume  (cm^/g):  0.99600 
Density  (g/cm^):  1.58000 
Input  Prop,  (kg):  1.00000 
Target  LD  (g/cm^);  0.95000 


RESULTS 


Velocity  (m/s):  1.709.2450 

Burnout  (cm):  157.9224 

Cham.  Vol.  (L):  10.50000 

Prop.  Mass  (kg):  9.97500 

Prop.  LD  (g/cm^):  0.95000 

Expansion  Ratio  (-):  6.1 1631 

Muzzle  Exit  Press  (MPa):  98.47 

Mass  Charge/Mass  Proj  (c/m  ratio):  1.10833 

Piezometric  Efficiency:  0.426 


Energy  Information 

Prop.  Energy  (MJ):  50.53999 
Proj.  KE  (MJ):  13.14683  (26.0%) 
GasKE(MJ):  4.85702  (9.6%) 

Gas  Internal  (MJ):  32.53614(64.4%) 


This  option  is  similar  to  Option  4  except  that  the  chamber  volume  is  taken  to  be  the  minimal  volume 
to  accommodate  the  charge  mass.  The  chamber  volume  is  computed  as  described  in  Option  3.  All 
comments  concerning  Option  4  apply  in  this  option  with  the  exception  of  the  chamber  volume  in  line  7. 
For  Option  5.  the  chamber  volume,  entered  on  line  7.  represents  the  maximum  acceptable  chamber  volume 
and  is  utilized  in  determining  the  maximum  charge  mass  considered  in  the  search. 
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Sample  input  and  ouqjut  for  Option  S  are  provided  in  Tables  14  and  IS. 


Table  14.  CONPRESS  Input  Deck  for  Options  5  and  6 


cpv21.ot5 

Test  of  the  output  routines  and  option  5 
5 

575. 

1.9. 

1.225  .996  1140.  1.58 
25.  120.  475. 

1500  .95 
7 


Table  IS.  CONPRESS  Output  for  Input  Deck  of  TaUe  14 


r 

INPUT  PARAMETERS 

i  Input  File:  cpv21.opS 

1  Run  Option:  5 

Output  File:  cpv21.ot5 

1  Gun  Parameters 

ProDCllani  Parameters 

1  Cham.  Vol.  (L):  25.0000 

Impetus  (J/g):  1,140.00 

1  Tube  Diam.  (mm):  120.0000 

Gamma  (-):  1.22500 

Proj.  Trav.  (cm):  475.0000 

Covolumc  (cmVg):  0.99600 

Proj.  Mass  (kg):  9.00000 

Density  (g/cm^):  1.58000 

Max  Press.  (MPa):  575.0000 

Input  Prop.  Ocg):  1.00000 

Target  VcL  (m/s):  liOO.0000 

Target  LD  (g/cm^:  0.95000 

RESULTS 

Ballistic  Valuc.s 

EnerRY  Information 

1  Velocity  (m/s):  1.773.6740 

Prop.  Energy  (MJ):  80.85132 

1  Burnout  (cm):  264.2473 

Proj.  KE(MJ):  14.15664  (17.5%) 

1  Cham.  Vol.  (L):  16.79735 

Gas  KE  (MJ):  8.36683  (TO.3%) 

1  Prop.  Mass  (kg):  15.95750 

1  Prop.  LD  (g/cm^):  0.95000 

1  Expansion  Ratio  (•)*  4.19820 

1  Muzzle  Exit  Press  (MPa):  151.01 

Gas  Ifticmal  (MJ):  58.32785  (72.1%) 

j  Mass  Charge/Mass  Proj  (cAn  ratio): 

1.77306 

1  Piezomctric  Efficiency:  0.458 
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To  validate  the  search  algorithm,  partial  results  from  a  parametric  run  of  CONHIESS  (Option  6)  in 
which  the  charge  mass  and  chamber  volume  were  varied  are  given  in  Table  16;  the  input  deck  is  shown 
in  Table  17.  As  can  be  seen  from  a  comparison  of  Tables  15  and  16,  the  correct  maximum  velocity  was 
obtained  utilizing  Option  5. 

TEST  CASE  6:  Multiple  Runs 

This  test  case  is  provide  to  illustrate  the  structure  of  the  input  deck  for  multiple  runs  of  COMPRESS. 
An  input  deck  for  two  runs  is  given  in  Table  18  with  the  corresponding  output  in  Table  19  (compare  to 
Tables  6  and  9). 


Table  16.  Partial  COMPRESS  Results  of  a  Parametric  Run  to  Determine  Optimal  Velocity 


Maximum  Pressure 
(MPa) 

Prop.  Mass 
(kg) . 

Cham.  Vol. 
(liters) 

Load.  Den. 
(g/cm^) 

Velocity 

(m/s) 

Burnout 

(cm) 

575.0000 

15.8000 

16.7000 

0.9461 

1,773.1320 

260.8842 

575.0000 

15.8000 

16.8000 

0.9405 

1.772.3870 

260.1642 

575.0000 

15.8000 

16.9000 

0.9349 

1,771.6400 

259.4406 

575.0000 

15.8500 

16,7000 

0.9491 

1.773.5420 

262.1747 

575.0000 

15.8500 

16.8000 

0.9435 

1.772.8010 

2f\4529 

575.0000 

15.8500 

16.9000 

0.9379 

1.772.0590 

260.7311 

575.0000 

15.9000 

16.8000 

0.9464 

1.773.2050 

262.7438 

1  575.0000 

15.9000 

16.9000 

0.9408 

1,772.4660 

262.0220 

1  575.0000 

15.9500 

16.8000 

0.9494 

1,773.5970 

264.0350 

I  575.0000 

15.9500 

16.9000 

0.9438 

1,772.8630 

263.3132 

1  575.0000 

16.0000 

16.9000 

0.9467 

1.772.2490 

264.6047 
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Table  17.  Input  Deck  That  Produced  the  Results  Shown  in  Table  16 


valopS.out 

Test  of  the  output  routines  and  Option  6 
6 

575. 

15.9575  9. 

1.225  .996  1140.  1.58 

10.5  120.  475. 

1500  .95 
7 


Table  18.  Example  of  COMPRESS  Input  Deck  for  Multiple  Runs 


multi.out 

This  tests  option  1  for  latest  version  of  COMPRESS 
1 

575. 

9.  9. 

1.225  .996  1140.  1.58 

10.5  120.  475. 

1200.  .95 

This  will  validate  option  2  of  latest  version  of  COMPRESS 
2 

575. 

1.9. 

1.225  .996  1140.  1.58 

10.5  120.  475. 

1549.07  .95 
7 
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Table  19.  COMPRESS  Results  for  Input  Deck  Shown  in  Table  18 


INPUT  PARAMETERS 

Ii^utFile:  multi.in 

Run  Option:  1 

Output  File:  multi.out 

Gun  Parameters 

Prooellant  Parameters 

Qtam.  Vol.  (L):  10.5000 

Tube  Diam.  (mm);  120.0000 

Proj.  Trav.  (cm):  475.0000 

Proj.  Mass  (kg):  9.00000 

Max  Press.  (I^):  575.0000 

Target  Vel.  (m/s):  1^00.0000 

Impetus  (J/g):  1,140.(K) 

Gamma  (-):  1,22500 

Covolume  (cnivg):  0.99600 

Density  (g/cm^):  1.58000 

Input  Prop,  (kg):  9.00000 

Target  LO  (g/cm^;  0.95000 

RESULTS 

Ballistic  Values 

Energy  Information 

Velocity  (m/s):  1,667.5470 

Bum  Out  (cm);  133.8050 

Cham.  Vol.  (L):  10.50000 

Prop.  Mass  (kg):  9.00000 

Prop.  LD  (g/cm^):  0.85714 
Expansion  Ratio  (-):  6.11631 
Muzzle  Exit  Press  (MPa):  88.31 
Mass  Charge/Mass  Proj  (c/m  ratio); 
Piczometric  Efficiency:  0.405 

Prop,  Energy  (IvU):  45.59999 

Proj.  KE  (MJ):  12.51320  (27.4%) 
GasKE(MJ):  4.17107(9.1%) 

Gas  Internal  (MJ);  28.91572  (63.4%) 

1.00000 

INPUT  PARAMETERS 

Input  File;  multi.in 

Run  Option;  2 

Output  File:  multi.out 

Gun  Parameters 

ProDcllant  Parameters 

Cham.  Vol.  (L):  10.5000 

Tube  Diam.  (mm);  120.0000 

Proj.  Trav.  (cm);  475.0000 

Proj.  Mass  (kg):  9.00000 

Max  Press.  (MPa):  575.0000 

Target  Vcl,  (m/s):  1349.0700 

impetus  (J/g);  1 . 1 40.00 

Gamma  (-):  1^2500 

Covolumc  (cni7g):  0.99600 

Density  (g/cm^);  1.58000 

Input  Prop,  (kg):^  1.00000 

Target  LD  (g/cm^);  0.95000 

RESULTS 

Balli.stic  Values 

Energy  Infomiation 

Velocity  (tn/s):  1.549.0700 

Bum  Out  (cm):  85.0453 

Cham.  Vol.  (L):  10.50000 

Prop.  Mass  (kg);  7.00109 

Prop.  LD  (g/cm^):  0.66677 
Ex{^sion  Ratio  (-);  6.1 1631 
Muzzle  Exit  Press  (MPa):  68.27 
Mass  Charge/Mass  Proj  (c/m  ratio): 
Piczometric  Efficiency:  0.350 

Prop.  Errergv  (MJ);  35.47216 

Proj  KE  (MJ):  10.79828  (30.4%) 
OasKE(MJ):  2.79999  (7.9%) 

Gas  Irttcmal(MJ):  21.87390(61.7%) 

0.77790 
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4.  CONCLUSIONS 


CONPRESS  is  a  constant  breech  pressure  (CBP)  interior  ballistic  code  which  not  only  perfonns  the 
basic  CBP  calculation  but  also  incoiporates  optimization  and  parametric  variation  algorithms.  In  this 
report  the  derivation  of  the  basic  CBP  equations  is  presented  together  with  a  user’s  manual  on  the  use  of 
the  program.  Validation  of  each  option  is  also  provided.  It  is  hoped  that  CONFUESS  will  prove  to  be 
a  useful  tool  in  evaluating  the  potential  of  new  propellant  formulations  and  novel  propulsion  concepts. 
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APPENDIX  A: 

COMPRESS  SOURCE  CODE  USTING: 
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PROGRAM  CONPRESS 

Q******************************************************************** 

C  Version  2,  December  1992 
C 

•  C  Author:  William  Oberle 
C  Advanced  Propulsion  Technology  Branch 

C  Weapons  Technology  Directorate 

'  C  U.S.  Army  Research  Laboratory 

C  Aberdeen  Proving  Ground,  Maryland,  21005“5066 

C  (410)  278  -  6200 

C 

C  This  program  will  compute  performance  for  a  constant  pressure  gun. 

C  It  has  several  options: 

C  1.  Determine  the  velocity,  KE  of  projectile  and  location  of  burn  out 
C  given  the  propellant  charge  mass  and  gun  geometry  including 
C  chamber  volume . 

C  2.  Determine  the  charge  mass  given  the  desired  velocity  and  chamber  vol. 

C  3.  Determine  the  charge  mass  and  chamber  volume  given  a  desired  velocity 
C  4 .  Determine  the  optimal  velocity  give  the  chamber  vol . ;  output  is  charge 
C  mass,  optimal  velocity  and  location  of  burn  out. 

C  5.  Determine  the  optimal  velocity  with  charge  mass  as  the  parameter, 

C  the  chamber  volume  is  taken  to  be  the  minimal  volume  needed  to  hold 
C  the  propellant  (subject  to  a  specified  loading  density) ;  output  is 
C  velocity,  chamber  volume,  location  of  burn  out  and  charge  mass. 

C  6.  Allows  parametric  runs  with  variations  in  all  input  parameters. 

C  The  calculation  can  be  performed  for  any  of  the  5  computational 

C  options . 

C  7.  Flag  for  termination  of  input  in  multiple  runs. 

QAtAAAlk*******'*******************************'************************** 

COMMON  /XLIST/  PC,MC,F,MP,GAMMA,COV,VC, A,TRAV,KINP, VEL,XB1,RH0, 

1  CVF,XLD,DVEL,IOPT,VCC,XMMC,D 

COMMON  /TIT/  TITLE, FILEIN,0UTF1LE 
CHARACTER*20  OUTFILE, FILEIN,TITLE*60 
REAL  MP,MC,KINP 
EXTERNAL  ALL,FUNC 
TOL-. 000001 

Qii*****ii***ii**1>**1iititii****1,i,*iiit*ii*ii1fk*****ii*1t1H(*1i*it***ii*****1t**1iii1i1i1tit*it1i 

Q************  First  the  input  file  name  is  determined  ****************** 

805  WRITE(*,800) 

800  FORMAT (//, 20X,  ••  CONPRESS,  Version  2  -  December  1992',///, 

IlOX,'  Enter  the  name  of  the  input  file.', 

2/,10X,'  (Enter  20  to  change  tolerance  used  in  search  routines.]',/) 
READ  ('*,801)  FILEIN 

801  FORMAT  (A20) 

Q**  **********  til***  ****■»*■»  **11* 

Q******************  Feature  to  allow  changing  tolerance  *************** 
C******************  in  option  4  fi  5  searches  ************************** 

Q*  *****************************************  ******i,**imt**************** 

IF  (FILEIN  .EQ.  '20')  THEN 
WRITE(*,851)TOL 

851  FORMAT (//,'  The  current  tolerance  for  searches  in  options' 

1, /,'  4  and  5  is:  ',F10.8,'.  Do  you  wish  to  enter  a  new  value' 

2, /,'  (Yes  -  1,  No  -  2)',/) 

READ(*, *) ITOL 

IF  (ITOL  .EQ.  1)  THEN 

WRITE (*,*) 'Enter  new  tolerance.' 
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READ(*,*)TOL 
END  IF 
GOTO  805 
END  IF 

Q************  ******  ******************************  ***★*****•*★****■*****★★ 

C****************  Test  that  file  exists  ******************************* 

OPEN {17 , FILE=FILEIN, STATUS=' OLD' , ERR=802) 

REWIND (17) 

GOTO  804 

802  WRITE (*, 803) FILEIN 

803  FORMAT(/,'  THE  SPECIFIED  INPUT  FILE',A20,'  DOES  NOT  EXIST!!', 

1/,'  PLEASE  CHECK  THE  FILE  NAME  A  ENTER  THE  CORRECT  FILE  NAME.', 

2///) 

PAUSE 
GOTO  805 

C********  Program  assumes  that  all  output  from  a  single  input  ********* 
C**************  deck  will  be  written  to  the  same  output  file  ********** 

804  READ(17,801)OUTFILE 
OPEN (UNIT=18, FILE-OUTFILE) 

C*******************  The  input  is  read  ******************************** 

Q'kicitir'k-ticit^iciticic'k'kir'kititiririr-kiric-k-kic-k'k'k'kicitic'k'k'kirir-kicirir'kifie'k'k-k'k'kic-kic'kiritic'k'k'k'k'k'k'kit'k'k'k 

1200  CALL  INPUT(IOPT,PC,MC,MP,GAMMA,COV,F,VC,A,TRAV,RHO,DVEL,CVF,VCC 

1,XMMC,D,XLD) 

Qieiriiitit'ki(iri(iti(^if  *  it*  icii  "kit  *  it  it  ic  it  it  if  if  irific  it  it  i(  it  ir  it  it  it  it  ir  it  it  irilic  it  iric  it 'k  if  it  it  it  it -kit -kirifir  "kit  ie  it 'k*  it  it 

C***************  The  correct  option  is  called  ************************* 


IF 

(lOPT 

.EQ. 

1) 

GOTO 

100 

IF 

(lOPT 

.EQ. 

2) 

GOTO 

200 

IF 

(lOPT 

.EQ. 

3) 

GOTO 

300 

IF 

(lOPT 

.EQ. 

4) 

GOTO 

400 

IF 

(lOPT 

.EQ. 

5) 

GOTO 

500 

IF 

(lOPT 

.EQ. 

6) 

GOTO 

600 

IF 

(lOPT 

•  EQ. 

7) 

GOTO 

1000 

C********  The  energy  of  the  projectile  is  computed  ************>•******* 

100  CALL  ENERGY(PC,MC,F,MP,GAMMA,COV,VC,A,TRAV,KINP) 

C********  The  velocity  and  location  of  burnout  is  computed  ************ 
CALL  XB ( P C , MC , F , MP , GAMMA, CO V , VC , A , XBl ) 

CALL  VELO{KINP,MP,VEL) 

Qiticititiiit*itiHiiHiit'kiiitiiitititHHHHititltii*^ititit'itiiitit1tititiritit1tifk>ititiiititiiitir1ti(1titititifHifiticitiii(it 

The  output  for  option  one  is  written 

CALL  XOUT 
GOTO  1200 


End  of  option  1 . 
Option  2 


AMC-VC*BH0‘ .1 
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BMC-VC* (1 .0-CVF) *RHO 

247  CALL  ENERGY (PC, AMC,F,MP, GAMMA, COV, VC, A, TRAV,KINP) 

AVEL»(2.0*KINP/MP) **.5-DVEL 

CALL  ENERGY (PC, BMC, F,MP, GAMMA, COV, VC, A, TRAV,KINP) 

BVEL* (2.0*KINP/MP)**. 5-DVEL 
IF( (AVEL*BVEL) .GT.0.0)  THEN 

WRITE (*,*) 'Bisection  Method  will  fail.' 

WRITE (*,*)' Enter  the  interval  to  be  searched.  OR' 
WRITE (*,*} 'Type  Control  C  to  terminate  the  program.' 
WRITE (*,*) 

WRITE (*,*)' Enter  the  endpoints  on  the  same  line.' 
WRITE (*,*) 

READ (*,*)AMC, BMC 
GOTO  247 

END  IF 

IF  (ABS(AVEL) .LE.  0.01)  THEN 
MC-AMC 

CALL  XB  (PC ,  MC ,  F ,  MP ,  GAMMA,  COV,  VC ,  A,  Xil ) 

CALL  VELO(KINP,MP,VEL) 

GOTO  99 

END  IF 

IF  (ABS(BVEL) .LE.  0.01)  THEN 
MC-BMC 

CALL  XB(PC,MC,F,MP,GAMMA,COV,VC,  A,XB1) 

CALL  VELO(KINP,MP,VEL) 

GOTO  99 

END  IF 

44  MC-.5*  (AMC+BMC) 

CALL  ENERGY (PC,  MC, F,MP, GAMMA, COV, VC, A, TRAV,  KINP) 

CVEL- ( 2 . 0  *KINP/MP ) *  * , 5-DVEL 

IF  (ABS(CVEL) .LE.  0.01)  THEN 

CALL  XB(PC,MC,F,MP,GAMMA,C0V,VC,A,XB1) 

CALL  VELO (KINP, MP,VEL)  . 

GOTO  99 
END  IF 

IF  ((AVEL*CVEL) .LE.0.0)  THEN 
BMC-MC 

ELSE 


AMC-MC 

END  IF 
GOTO  44 

99  CALL  XOUT 

GOTO  1200 

'ff  if  'k  'k  'k  it  "k  it  'k  'k  it  "k  ^  ii  -k  -k  "k  -k  it  *  it  it  it  it  ir  it  ic  i.  -tc  it  ir  it  it  it  ^  it  it  it  it  it  it  it  it  it  it  It  it  -k  it  -it  it  it  it  it  it 

C*********************  option  2  ’■'«f*'>«****'»*'»*********i********'»*'* 

C*  ***************************************  It* 

C************************  Option  3  *********************************** 

C 

C*******  The  initial  charge  mass  is  the  minimum  mass  expected.  ******** 
C*******  A  search  will  be  pertormed  to  determine  the  endpoints  ******** 
C*******  of  the  interval  for  vropellant  masses  starting  with  *♦**★♦★* 
C*******  the  initial  charge  Mass.  ♦*♦**★** 

Q**  **********************^r->^  ******  it****it******ittt  ******  it*  Itiiiiit*ii*1iviiiniitii1i1i 

300  NTEST-1 

AMC-MC 
VC-MC/RHO 
VC-VC/(1. 0-CVF) 

CALL  ENERGY (PC, AMC, F, MP, GWiMA, COV, VC, A, TRAV, KINP) 
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AVEL- ( 2 . 0  *KINP /MP ) *  * . 5-DVEL 
BMC-AMC 

111  BMC«BMC*1.2 
VC-BMC/RHO 
VC-VC/(1.0-CVF) 

CALL  ENERGY (PC, BMC, F, MP , GAMMA, COV, VC, A, TRAV, KINP) 

BVEL- (2 . 0*KINP/MP) ** . 5-DVEL 
IF ( (AVEL*BVEL) . GT . 0 . 0 )  THEN 
NTEST-NTEST+1 
IF  (NTEST  .GT.  500)  THEN 
WRITE (*,112) 

112  FORMAT ('  A  suitable  right  hand  endpoint  has  not  been  found', 
1/,'  after  500  iterations.  This  corresponds  to  multiplying  ', 

2/,'  the  input  value  of  the  propellant  mass  by  3E+36.  It  ', 

3/,'  appears  that  the  desired  velocity  cannot  be  achieved  in', 

4/, 'the  given  projectile  travel.  Run  option  5  to  determine  if' 
5/,'  the  desired  velocity  can  be  achieved.',//) 

PAUSE 
GOTO  1000 
END  IF 
GOTO  111 
END  IF 

Q***********************************************************i»********** 

C*******  The  endpoints  of  the  interval  have  been  determined  *********** 

Qicir  icir  t  Tkikicitic  ic  it  it  irititir  it  itiririric  ir  it  irir  it  it  it  icitiriririrititititikititiririririrititTk’kiritWitititiriritiriririrititiritii 

IF  (ABS(AVEL) .LE.  0.01)  THEN 
MC-AMC 
GOTO  199 

ENDIF 

IF  (ABS(BVEL) .LE.  0.01)  THEN 
MC-BMC 
GOTO  199 

ENDIF 

144  MC-.5* (AMC+BMC) 

VC-MC/RHO 

VC"VC/(1.0-CVF) 

CALL  ENERGY (PC, MC, F, MP, GAMMA, COV, VC, A, TRAV, KINP ) 

CVEL- (2 . 0*KINP/MP) ** , 5-DVEL 
IF  (ABS(CVEL) .LE.  0.01)  GOTO  199 
IF  ((AVEL*CVEL) .LE.O.O)  THEN 
BMC-MC 
GOTO  155 

ENDIF 

AMC-MC 

155  GOTO  144 

199  VC-MC/RHO 

VC-VC/ (1.0-CVF) 

CALL  XB (PC, MC, F,MP, GAMMA, COV, VC, A, XBl ) 

CALL  VELO(KINP,MP,VEL) 

CALL  XOUT 
GOTO  1200 

^****<**«****«************«****«**«**«**«*«*********«****««*****«««***1r* 

C********************  End  of  option  3  ********************************* 
(2*  *********************************************  **********************  * 
Q*** ***********************************************************  ******** 

Q*************************  Option  4  *********************************** 
0  ********************************************************************* 
C*******  The  values  to  bracket  the  charge  mass  are  determined.  ******* 
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C*******  The  max  value  for  charge  mass  is  the  max  that  can  be 
Q*******  held  in  the  chamber  accounting  for  loading  density 

400  AX-0.0 

c************************************************************** 

C********************  AX  is  left  hand  endpoint  **************** 
C***********  cx  is  maximum  amount  which  can  be  held  in  chamber 
C***********  Will  determine  the  velocity  at  endpoint  CX  since 
Qic **********  the  search  routines  sometimes  miss  the  endpoint 

Q************************ ************************************** 

CX-VC*(1.0-CVF) *RHO 

CALL  ENERGY (PC, CX, F, MP, GAMMA, COV, VC, A, TRAV, KINP ) 

CALL  XB (PC , CX , F , MP , GAMMA, COV, VC, A, XBl ) 

CALL  VELO(KINP,MP,VEL) 

VELT-VEL 

KINPT-KINP 

XBIT-XBI 

MCT-CX 

BX-CX/2.0 

VEL-BRENT (AX, BX, CX, ALL, TOL, XMIN) 

VEL— 1.0*VEL 

IF  (VEL  .LT.  VELT)  THEN 

VEL-VELT 

MC-MCT 

XBI-XBIT 

MC-MCT 

ELSE 

MC-XMIN 

ENDIF 

CALL  XOUT 

GOTO  1200 

Q**********************  £nd  of  option  4  *********************** 

Q*** ***************************************************  ******** 

Q**********************  Option  5  ****************************** 

Q******««*****  <j>he  input  chamber  volume  is  assumed  to  be  ****** 
0*************  the  largest  acceptable  volume  for  the  run  ****** 

Q******  **************************************************  ****** 

500  AX-0.0 

CX-VC* (1 .0-CVF) *RHO 

Q************************************************************** 

Q***********  Q}(  maximum  amount  which  can  be  held  in  chamber 
C***********  Will  determine  the  velocity  at  endpoint  CX  since 
Q***********  the  search  routines  sometimes  miss  the  endpoint 

£••****•*•**•*••*•«**•••*•«*•*•*•***«****•*•*•*****•********•*« 

CALL  ENERGY (PC, CX, F,MP, GAMMA, COV, VC, A, TRAV, KINP) 

CALL  XB (PC,  CX, F, MP , GAMMA, COV, VC, A, XBl ) 

CALL  VELO (KINP, MP, VEL) 

VELT-VEL 

KINPT-KINP 

XBlT-XBl 

MCT-CX 

BX-CX/2.0 

VEL-BRENT (AX, BX, CX, FUNC, TOL,  XMIN) 

VEL— 1 .0*  VEL 

IF  (VEL  .LT.  VELT)  THEN 

VEL-VELT 


★★★★★★★ 

******** 

******** 

******** 

******** 

******** 

******** 

******** 


******** 
******  ** 
******** 
******** 
******** 
******** 
******** 
******** 


******** 

******** 

******** 
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MC-MCT 

XBI-XBIT 

MC-MCT 

ELSE 

MC-XMIN 

ENDIF 

VOL-MC/RHO 
VOL-VOL/ (1.0-CVF) 

CALL  XOUT 
GOTO  1200 

Q *★★★★******★**★★*★**★**★*★★*★******★**★★★*******★★★*★*★*★★★★★★★★★*★*** 

C**********************  End  of  option  5  ******************************* 
C***********************  OPTION  6  ************************************* 

600  WRITE (*,601) 

601  FORMAT(//,'  This  option  will  allow  parametric  runs  of  CONPRESS:' 

1, /,'  Starting,  ending  and  incremental  values  for  each  variable' 

2, /,'  will  be  requested.  If  the  variable  is  not  to  be  included' 

3, /,'  in  the  parametric  search,  the  input  data  value  will  be  used.' 

4, /,'  Output  will  be  a  file  of  values  from  the  parametric  run.') 

PAUSE 

^  if 'k  it -k -k  ir  it  ii  it  if  it  it  ic -k  it  it  -k  it  it  it  it  "k  it  it  it  it  it -it  ic  it  ^r  it  it  "k  it  it  it  it  it  ir -k it  it  it  "k  it  it  it  it  if  it  it  ir  it  i(  it  it  it  "k  it  ir  if  it  it  it  i: 

C****************  CALCULATION  OPTION  IS  DETERMIND  ★★★★*★**★***★★******* 

WRITE (*,*)' Enter  the  calculation  option  (1  -  5).' 

READ(*, *) IOP6 
WRITE (18, 607) IOP6 

607  FORMAT ('l',20X, 'CONPRESS  VERSION  2,  December  1992', 

1//,15X, '  A  parametric  search  is  being  performed.', 

2/,15X,'  The  computational  option  is;  ',11,//) 

Q******-k*******1i***-klfk***1i-tit*iHi1t**1iiHi*1i*itit****iHt**1t*1)i>t**************** 

Q********^*********  Output  file  headers  are  written  ******************* 

0****«****«**********«***lir**************<r*«****'»*************«********* 

WRITE(18,699) 

699  FORMAT(/,'  Max  Press  Prop  Mass  Proj  Mass  Gamma  Covol  Imp 

letus  Density  Cham  Vol  Tube  Diam  Travel  Load  Den  Velocity 
2  Proj  KE  Burn  Out  ') 

WRITE (18, 795) 

795  FORMAT(  '  (MPa)  (kg)  (kg)  (-)  (cc/g)  (J 

1/g)  (g/cc)  (Liters)  (mm)  (cm)  (g/cc)  (m/s) 

2  (MJ)  (cm)') 

O****'^*****'***************^**'**************************'**************** 

WRITE (*,*) 'Vary  chamber  pressure?  (Yes  -  1,  No  «  2)' 

READ(*,*>ITEST 

IF  (ITEST  .EQ.  1)  THEN 

WRITE(*,602) 

602  FORMAT (/,'  Enter  the  starting  value,  ending  value  and  increment' 

1, /,'  all  on  the  same  line  separated  by  a  space.' 

2, /,'  Units:  Pressure  in  MPa' 

3, /,'  Charge  &  Projectile  Mass  in  kg' 

4, /,'  Covolume  in  cc/g' 

5, /,'  Impetus  in  J/g' 

6, /,'  Density  in  g/cc' 

7, /,'  Chamber  Volume  in  liters' 

8, /,'  Tvibe  Diameter  in  mm' 

9, /,'  Projectile  Travel  in  cm' 

*,/,'  Loading  Density  in  g/cc') 
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READC*,  *)SPC,EPC,DPC 

Q***********-***  Conversion  to  MKS  units  for  program  ******************** 
C***************  mKS  needed  for  loop,  1E7  ARE  TOO  BIG  ****************** 
C**********  Will  convert  to  cgs  at  the  loop  **************************** 
ELSE 

SP>PC/1.E7 
EPC“PC/1.E7 
DPC=PC/1.E7 
END  IF 

Q^*****ilf*Hf**w***^i*r***-************************************************** 


IF  (I0P6  .NE.  1)  THEN 
SMC-MC 
EMC-MC 
DMC-MC 
GOTO  1487 
END  IF 

WRITE {*,*)' Vary  charge  mass?  (Yes  -  1 

READ(*,*)ITEST 

IF  (ITEST  .EQ.  1)  THEN 

WRITE (*, 602) 

READ(*, *)SMC,EMC,DMC 

C**************  Conversion  to  grams  for  program 
SMC- SMC *1000. 


,  No  -  2)' 


■k*’k**i(***ifk***  ■',******** 


EMC“-EMC*1000. 

DMC-DMC*1000. 

ELSE 

SMC-MC 

EMC-MC 

DMC-MC 

ENDIF 


1487  WRITE {*,*) 'Vary  projectile  mass?  (Yes  -  1,  No  -  2)' 

READ(*,  *)  ITEST 

IF  (ITEST  .EQ.  1)  THEN 

WRITE{*,602) 

READ(*,  '‘)SMP,EMP,DMP 

C**************  Conversion  to  grams  for  program  *********************** 

SMP-SMP*1000. 

EMP-EMP*1000. 

DMP-DMP*1000, 

ELSE 

SMP-MP 

EMP-MP 

DMP-MP 

ENDIF 


WRITE (*,*)' Vary  gamma?  (Yes  -  1,  No  -  2)' 

READ(*, *) ITEST 

IF  (ITEST  .EQ.  1}  THEN 

WRITE(*,602) 

READ ( * , * ) SGAM, EGAM, DGAM 
ELSE 

SGAM- GAMMA 
EGAM-GAMMA 
DGAM-GAMMA 
ENDIF 

C*****************<**************************************************** 

WRITE {*,*) 'Vary  covolume?  (Yes  -  1,  No  -  2) ' 
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READ(*,*)ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ ( * , * ) SCOV, ECOV, DCOV 
ELSE 

SCOV^COV 
ECOV-COV 
DCOV=COV 
END  IF 

C***»******************  ************************************************* 

WRITE (*,*) 'Vary  propellant  impetus?  (Yes  -  1,  No  «  2)' 

READ (*,*) ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ(*, *)SF,EF,DF 

C**************  Conversion  to  cgs  units  for  program  ******************* 
SF-SF*l.E+07 
EF-EF*l.E+07 
DF-DF*1 .E+07 
ELSE 
SF=F 
EF-F 
DF-F 


END  IF 

Q****************************************************** **************** 

WRITE (*,*) 'Vary  propellant  density?  (Yes  -1,  No  -  2)' 

READ(*, *) ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ ( * , * ) SRHO, ERHO,  DRHO 

ELSE 

SRHO-RHO 

ERHO-RHO 

DRHO-RHO 

ENDIF 


IF  ((IOP6  .EQ.  3)  .OR.  (IOP6  .EQ.  5))  THEN 
SVC-VC 


EVC-VC 
DVC-VC 
GOTO  9842 
ENDIF 

WRITE (*,♦) 'Vary  chamber  volume?  (Yes  -  1,  No  -  2)' 
READ(*, *) ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ(*, *)SVC,EVC,DVC 


Q**************  Conversion  to 
SVC-SVC*1000. 
EVC-EVC*1000. 
DVC-DVC*1000. 

ELSE 

SVC-VC 

EVC-VC 


cc  for  program  •••****♦*••**••***•**•***• 


DVC-VC 

ENDIF 


9842  WRITE (*,*) 'Vary  tube  diameter?  (Yes  -  1,  No  -  2)' 
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READ(*,*)ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ(*,*)SD,ED,DD 

C**************  Conversion  to  cm  for  program  ************************** 
SD-SD/10. 

ED-ED/10. 

DD-DD/10 . 

ELSE 

SD-D 

ED-D 

DD-D 

ENDIF 

Q**************r*****'(t**'*r**i»**'**'i(*************************************** 

WRITE (*,*)' Vary  projectile  travel?  (Yes  -  1,  No  -  2)' 

READ C*,*) ITEST 
IF  (ITEST  .EQ.  1)  THEN 
WRITE (*,602) 

READ ( * , * ) STRAV, ETRAV, DTRAV 
ELSE 

STRAV-TRAV 

ETRAV-TRAV 

DTRAV-TRAV 

ENDIF 

Q********«************i^*«****«****** *************** ******************** 

IF  (IOP6  .EQ.  1)  THEN 

SXLD-XLD 

EXLD-XLD 

DXLD-XLD 

GOTO  3876 

ENDIF 

WRITE (*,*) 'Vary  loading  density?  (Yea  -  1,  No  -  2)' 

READ (*,*) ITEST 
IF  (ITEST  .EQ.  i)  THEN 
WRITE (*,602) 

READ ( * , * ) SXLD , EXLD , OXLD 
ELSE 

SXLD-XLD 

EXLD-XLD 

DXLD-XLD 

ENDIF 

(1^**  ************************************************************  *•***«*• 

3876  DO  650  PC1-SPC,EPC,DPC 
PC«PC1*1.E7 

DO  650  MC1-S«C,EMC,DHC 
DO  650  MP-SMP,EMP,DMP 
DO  650  GAMMA-5GAM,EGAM,DGAM 
DO  650  COV-SCOV,ECOV,DCOV 
DO  650  F-SF,EF,DF 
DO  650  RHO-SRHO,ERHO,DRHO 
DO  650  VCl-SVC, EVC, DVC 
DO  650  D-SD,ED,DD 
A-3.i415926S359* (D/2.0) **2.0 
DO  650  TRAV-STRAV, ETRAV, DTRAV 
DO  650  XLD-SXLD,EXLD,DXLD 
CVF-1 .-(XLD/RHO) 

^•••******************^.’<t****»**************««*****«******************** 

C*****************  OPTICW  1  CALCOLATION  ******************************* 
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Qir  ir  it  it  it  ic  it  "k  ic  it  it  i:  it  it  ic  it  it  it  it  it  it  ir  ic  if  it  ie  if -k  ic -kit  it  it  it  it  it  ie  icit  "k  "ic  it  it  it  if  it  ir  it  it  it  it  it  "k  it  it  it  it  it  it  it  it 'k -kit  it  it  it  it  it 

IF  (I0P6  .EQ.  1)  THEN 

MC-MCl 

VC«VC1 

CALL  ENERGY (PC, MC, F , MP , GAMMA, COV, VC, A, TRAV, KINP ) 

CALL  XB (PC , MC , F , MP , GAMMA, COV, VC , A , XBl ) 

CALL  VELO(KINP.-MP,VEL) 

GOTO  631 
END  IF 

(;***★***■*****★********•******★*******■*****★**********★**■***★**★**★**★★** 
C*****************  OPTION  2  CALCULATION  ******************************* 

Qit  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  icic  it  it  it  it  It  it  icit  it  it  it  it  it  it  it  icit  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it ’k  it  it  irit  it  it  it  it  it  it  it 'i  it 

IF  (IOP6  .EQ.  2)  THEN 
VC-VCl 

AMC-VC*RHO* . 1 
BMC-VC* (1.0-CVF) *RHO 

2247  CALL  ENERGY (PC, AMC, F,MP, GAMMA, COV, VC, A, TRAV, KINP) 

AVEL- (2 .0*KINP/MP) ♦* . 5-DVEL 

CALL  ENERGY (PC, BMC, F,MP, GAMMA, COV, VC, A, TRAV, KINP) 

BVEL- (2 .0*KINP/MP) **. 5-DVEL 
IF( (AVEL*BVEL) .GT.0.0)  THEN 

WRITE (*,*) 'Bisection  Method  will  fail.' 

WRITE (*,*) 'Enter  the  interval  to  be  searched.  OR' 
WRITE (*,*) 'Type  Control  C  to  terminate  the  program. 
WRITE (*,*) 

WRITE (*,*) 'Enter  the  endpoints  on  the  same  line.' 
WRITE  (*, *) 

READ  (*,*)  fUMC,  BMC 
GOTO  2247 

END  IF 

IF  (ABS(AVEL) .LE.  0.01)  THEN 
MC-AMC 

CALL  X0 (PC, MC, r,MP, GAMMA, COV, VC, A, XBl) 

CALL  VEL0(KINP,MP,VEL) 

GOTO  631 

ENOIF 

IF  (ABSOVEL)  ,LE-  0.01)  THEN 
HC-BMC 

CALL  XB  ( PC,  HC ,  F,  MP ,  GAt-BiA,  COV,  VC ,  A ,  XBl ) 

CALL  VELO(KINP,MP,VSL) 

GOTO  631 

ENDIF 

244  MC-.5*  (A^iC*BMC) 

CALL  ENERGY  (PC, MC,  F,MP,  GAEm,COV,  VC,  A,TRAV,  KINP) 

CVEL-(2.0*KIKP/MP) *«.5~DV£L 

IF  (ABSfCVEL) .LE.  0.01)  THEN 

CALL  XB  (PC,  KC ,  F ,  MP ,  GAi'-tMA,  COV,  VC ,  A,  XQl ) 

CALL  VELO{KlNP,M?,VEL) 

GOTO  631 
ENDIF 

IF  ( (AVEL*CVEL) .LE.0.0)  THEN 
BKC-MC 

ELSE 


AMC*MC 


C 


ENDIF 
GOTO  244 
ENDIF 
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c*****************  OPTION  3  CALCULATION  ***********************-******** 

Q************ *****<■  **************** *********************  +  ************** 

IF  (IOP6  .EQ.  3)  THEN 

MC=MC1 

NTRST=1 

AMC==MC 

VC=MC/RHO 

VC=VC/ (1.0-CVF) 

CALL  ENERGY (PC, AMC, F,MP, GAMMA,COV, VC, A, TRAV, KINP) 

AVEL= (2 . 0*KINP/M:0  ** . 5-DVEL 
BMC=AMC 

3111  BMC-BMC*1,2 
VC=BMC/RHO 
VC=VC/ (1.0-CVF) 

CALL  ENERGY (PC, BMC.F,MP,GAMMA,COV, VC, A, TRAV, KINP) 

B\'EL=  (2.0  *KINP/MP )  *  * .  5-DVEL 
IF(  (A’/EL*BVEL)  .GT.0.0)  THEN 
NTEST=NTEST+1 
IF  (NTEST  .GT.  500)  THEN 
WRITE (*,3112) 

3112  FORMAT ('  A  suitable  right  hand  endpoint  has  not  been  found', 
1/,'  after  500  iterations.  This  corresponds  to  multiplying  ', 

2/,'  the  input  value  of  the  propellant  mass  uy  3E+36.  It  ', 

3/,'  appears  that  the  desired  velocity  cannot  be  achieved  in', 

4/, 'the  given  projectile  travel.  Run  option  5  to  determine  if' 

5/f '  the  desired  velocity  can  be  achieved.',//) 

PAUSE 
GOTO  ICOO 
END  IF 
GOTO  3111 
END  IF 

C**X********  *************.',•*****************  A*********-^***************** 

C*******  The  endpoints  of  the  interval  have  been  determined  *********** 
C************************«*********  *****,>*****************★************ 

IF  (ABS(AVEL) .LE.  O.Ol)  THEN 
MC^AMC 
GOTO  3199 

END  IF 

IF  (ABS(BVEL) .LE.  0.01)  THEN 
MC»BMC 
GOTO  3199 

ENDIF 

3144  MC-..5*  (AMC+BMC) 

VC-MC/RHO 
VC=VC/ (1.0-CVF) 

CALL  EN'  '••GY(PC,MC,F,MP,GAMMA,C0V,VC,A,TRAV,KINP) 

CVEL- (2 .C*KINP/MP) ** . 5-DVEL 
IF  (ABS(CVEL) .LE.  0.01)  GOTO  3159 
IF  ( (AVEL-CVEL) , LE . 0 . 0)  THEN 
BMC-MC 
GOTO  3155 

ENDIF 

AMC-MC 

3155  GOTO  3144 

3199  VC«MC/RHO 

VC^VC/ (1.0-CVF) 

CALL  XB (PC , MC , F , MP , GAMMA, COV, VC ■ A, XBl ) 

CALL  VELO(KINP,MP,VEL) 
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GOTO  631 
END  IF 

C*****************  OPTION  4  CALCULATION  ******************************* 

IF  (IOP6  .EQ.  4)  THEN 

VC=VC1 

AX“0 .0 

C********************  AX  is  left  hand  endpoint  ************************ 
C***********  cx  is  maximum  amount  which  can  be  held  in  chamber  ******** 
C***********  Kill  determine  the  velocity  at  endpoint  CX  since  ******** 
C***********  the  search  routines  sometimes  miss  the  endpoint  ******** 


CX=VC* (1.0-CVF) *RHO 

CALL  ENERG Y { P  C , CX , F , MP , GAMMA , CO V , VC , A , TRAV , KI NP ) 
CALL  XB (PC, CX, F, MP , GAMMA, COV, VC, A, XBl ) 

CALL  VELO{KINP,MP,VEL) 

VELT”VEL 

KINPT>=KINP 

XBIT-XBI 

MCT«=CX 

BX»CX/2.0 

VEL-^'BRENT  (AX,  BX,  CX,  ALL,  TOL,  XMIN) 

VEL“-1,0*VEI. 

IF  (VEL  .LT.  VELT)  THEN 

VEL=VELT 

MC-MCT 

XB1«XB1T 

MC-MCT 

ELSE 


MC-XMIN 
ENDIF 
GOTO  631 
ENDIF 

C.***A**. »****«*#*  OPTION  5  CALCULATION  ******************************* 


IF  (IOP6  .EQ.  5)  THEN 

VC-VCl 

AX-0.0 

CX-VC* (1.0-CVF) *RHO 


0i'********«*****«»*'«****4«****  *«•**«••«••*•••*«*•«*»•••«•******•***•**• 


C***********  OX  is  maximum  amount  which  can  be  held  in  chamber  ******** 
0***»*»***»«  Will  determine  the  velocity  at  endpoint  CX  since  ******** 
0*«***4***»*  the  search  routines  sometimes  miss  the  endpoint  ******** 

O************************************* *•«••*•*•*•«••*•**••••**•*«*••••* 


CALL  ENERGY ( PC , CX , F , MP , GAMMA , COV, VC , A , TRAV, Kl NP ) 
CALL  XB (PC, CX , F , MP , GAMMA, COV, VC, A, XBl ) 

CALL  VELO(KINP,MP,VEL) 

VELT- VEL 

KINPT-KINP 

XBlT-XBl 

MCT-CX 

BX-CX/2.0 

VEL-BRENT (AX, BX , CX, FUNC, TOL, XMIN) 

VEL— 1 .0*VEL 

IF  (VEL  .LT.  VELT)  THEN 


46 


o  n  n'n  oooooooooo 


VEL-VELT 

MC=MCT 

XB1=XB1T 

MC=MCT 

ELSE 

MC»XMIN 

END  IF 

VOL=MC/RHO 
VOL-VOL/ (1.0-CVF) 

GOTO  631 
END  IF 

Q******************  CONVERSION  FOR  PRINTOUT  *************************** 

631  Al-PC*lE-7 

A2-MC/1000. 

A3-t^P/1000. 

A4 -GAMMA 

A5-COV 

A6-F*1.E-7 

A7-RHO 

A8-VC/1000. 

A3-D*10. 

AlO-TRAV 

A11-A2/A8 

A12-VEL/100. 

A13-KINP*1.E-7*1.E-6 

A14-XB1 

C*********************  Values  written  to  file  ************************* 

WRITE  (18,  651)A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,A11,A12,A13,A14 
651  FORMAT  {14n0. 4) 

650  CONTINUE 

GOTO  1200 

Q****it*************'****'*****w*«**«***4r*«******************«************ 

C********************  end  of  option  6  ********************************* 

1000  close (UNIT-17) 

CLOSE (UNIT-18) 

END 

Q***«*'**A*«*****1»**«********«****«***«*«******«**«**i»*******4«**'******* 

C*******«**************«<r  Input  Subroutine  **************************** 

Q*«*****«**«**t'*****«*««««***i»**4*«*«**«*«***********  *-•*•««!»**••**««*•* 

SUBROUTINE  INPUT ( lOPT, PC, MC, MP, GAMMA, COV, F, VC, A, TRAV, RHO, DVEL, 

1  CVF, VCC,XMMC,D,XLD) 

***ft**********«*«*4t******«***«*«*******«*****«**«*'»«**i»**«i»*****««**** 

«*«******«*•«•*'«•*«**  Information  on  input  parameters  **************** 

**t*«**«*««**«***««****«**«***«**«**  «*••**«•*«••**•*••*  «’•*«•***•*«*'*•* 

INPUT  PARAMETERS  REQUIRED 

lOPT  —  Program  option  flag 

PC  —  Pressure  in  MPa 

MC,XMMC  —  Charge  mass  in  kg 

MP  —  Projectile  mass  in  kg 

GAMMA  —  Ratio  of  specific  heats 

COV  —  Propellant  covolume  in  cc/g 

F  —  Propellant  impetus  in  J/g 

RHO  —  Propellant  density  in  g/cc 

VC,VCC  —  Chamber  volume  without  propellant  in  L 
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D  —  Gun  tube  diameter  (ram) 

TRAV  —  Projectile  travel  in  cm 

DVEL  —  Desired  velocity  (m/s) 

XLD  —  Desired  propellant  loading  density  (g/cc) 

COMPUTED  VALUES 

A  —  Bore  cross-sectional  area  cm''2 

CVF  —  Chamber  void  fraction;  fraction  of  chamber 

which  cannot  contain  propellant 

iticiticifificiticicitifiri^i^i^iciciticiciciciciriritic-k-k'k'kik’k'k'kivif'kit'ic'k'kic'k'kicicic'k'kicit'k'k'kic'ic-k'kic'kic'kicir'k^ic'k 

COfiMON  /TIT/  TITLE,  FILEIN,OUTFILE 
CHARACTER*20  TITLE*60,FILEIN, OUTFILE 
REAL  MC,MP 

Q**  *************************************************************  ********** 

READ  (17, 800, END*=900)  TITLE 
800  FORMAT (A60) 

IF  (TITLE  .EQ.  '7')  GOTO  905 

READ (17, *, END-900 ) lOPT 

IF  (lOPT  .EQ.  7)  GOTO  905 

READ (17, *, END-900) PC 

READ (17,*, END-900 )MC,MP 

READ ( 17 , * , END-900 ) GAMMA, COV,  F, RHO 

READ (17, *, END-900) VC, D, TRAV 

READ ( 17 , * , END-900 ) DVEL, XLD 

C********************************************************************** 

C*******  Conversion  of  all  quantities  to  cgs  system  of  units  ********** 

PC-PC*1.0E+07 

MC-MC*1000.0 

XMMC-MC 

MP-MP*1000.0 

F-F*l .OE+07 

VC-VC*1000. 

VCC-VC 

D-D/10. 

A-3. 14159265359* (D/2.0) **2.0 
DVEL-DVEL*100. 

C****************  Compute  chamber  void  fraction  *********************** 

C********************************************************************** 

CVF-l.-(XLD/RHO) 

GOTO  901 

900  WRITE(18,902) 

902  FORIiAT  ( /  / , '  Unexpected  end  of  file  has  occurred.', 

1/,'  Program  Terminated') 

905  CLOSE (UNIT-17) 

CLOSE (UNIT-18) 

STOP 

901  RETURN 
END 

c********************************************************************** 
(;***#*#********<. **.*•  Subroutine  Energy  ****************************** 

c********************************************************************** 

SUBROUTINE  ENERGY (PC, MC, F, MP , GAMMA, COV, VC , A, TRAV,  KINP ) 

REAL  MC,MP,KINP 

V1F-VC-C0V*MC 

VMF-VC+A*TRAV-COV*MC 


CON1-1+MC/(2.0*MP) 

CON2-1+MC/(3.0*MP) 

C0N3-GAMMA-1 . 0 

VBF-MC*F*C0N1/ ( GAMMA* C0N2*PC) +C0N3*V1F/GAMMA 
C0N4- (VBF/VMF) **GAMMA 

KINP-PC/CON1*(GAMMA/CON3*VBF-V1F-VMF*CON4/CON3) 

RETURN 

END 

C********************************************************************** 

C*********************  Subroutine  XB  ********************************** 

SUBROUTINE  XB (PC, MC, F, MP , GAMMA, COV, VC, A, XBl ) 

REAL  MC,MP 

V1F-VC-C0V*MC 

CON1-1+MC/(2.0*MP) 

CON2-1+MC/(3.0*MP) 

CON3-GAMMA-1 . 0 

VBF-MC*F*CONl/ (GAMMA*CON2*PC) +C0N3*V1F/6AMMA 

VBFl=VBF+COV*MC-VC 

XBl-VBFl/A 

RETURN 

END 

^★★★★★★★★★★★★★★★★★vnt*************************************************** 

C*********************  Subroutine  VELO  ******************************** 

SUBROUTINE  'v/ELO  (KINP,MP,  VEL) 

REAL  MP,KINP 
VEL'“{2.0*KINP/MP)  **,5 
RETURN 
END 

C********  The  function  all  is  used  in  the  subroutine  BRENT  ************ 

Qi********************************************************************** 

FUNCTION  ALL { XX) 

COMMON  /XLIST/  PC, MC, F, MP, GAMMA, COV, VC, A, TRAV, KINP, VEL, XBl, RHO, 
1  CVF,XLD,DVEL, IOPT,VCC,XMMC,D 

REAL  MC,MP,KINP 
V1F-VC-C0V*XX 
VtiF-VC+A*TRAV-COV*XX 
CONl-l+XX/ (2.0*MP) 

CON2-1+XX/ {3.0*MP) 

CON3-GAMMA-1.0 

VBF-XX*F*CONl/ ( GAMMA *CON2* PC) +C0N3*V1F/GAMMA 
CON4- (VBF/VMF) **GAMMA 

KINP-PC/CONl* (GAMMA/CON3*VBF-V1F-VMF*CON4/CON3) 

VBFl-VBF+COV*XX-VC 

XBl-VBFl/A 

VEL-(2.0*KINP/MP) *».5 
ALL— 1,0*  VEL 
RETURN 
END 

Q*  *  *  ■»»■»*•»*■»*■»■»  H  *****  *i>  •  *  *  »  *  K  *  9  ****■»■»*  -t*  t  ***■»*«***■»*■»■*  »  t  m  *********  *****  * 

Q*********************  function  brent  *******«************************* 

Q**  ******************************************************************  ** 

function  brent ( ax, BX, CX, F, TOL, XMIN) 

EXTERNAL  F 

PARAMETER  ( ITMAX-100, C60LD- . 3819660, ZEPS-1 . OE-10) 

A-MIN(AX,CX) 
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B-MAX  (AX,  CX) 

V-BX 

W-V 

X-=V 

E-0. 

FX-F(X) 

FV^FX 

FW=FX 

DO  11  ITER-1,ITMAX 
XM»0.5* (A+B) 

TOLl=TOL*ABS (X) +ZEPS 
TOL2«2 . *TOLl 

IF(ABS(X-XM) .LE, (TOL2-.5*(B-A)) )  GOTO  3 
IF ( ABS (E) . GT . TOLl )  THEN 
R-(X-W)*(FX-FV) 

Q-(X-V)*(FX-FW) 

P-(X-V)*Q-(X-W)*R 

Q-2.*(Q-R) 

IF(Q,GT.O.)  P— P 
Q-ABS (Q) 

ETEMP-E 

E-D 

IF (ABS (P) .GE.ABS(.5*Q*ETEMP) .OR.P.LE.Q* (A-X) .OR. 

*  P.GE.Q*{B-X} )  GOTO  1 
D-P/Q 
U-X+D 

IF(U-A.LT.TOL2  .OR.  B-U.LT.TOL2)  D-SIGN (TOLl, XM-X) 
GOTO  2 
END  IF 

1  IF(X.GE.XM)  THEN 
E-A-X 

ELSE 
E-B-X 
END  IF 
D-CGOLD*E 

2  IF(ABS(D) .GE.TOLl)  THEN 
U-X+D 

ELSE 

U-X+SIGN(TOLl,D) 

END  IF 
FU-F(U) 

IF(FU.LE.FX)  THEN 

IF(U.GE.X)  THEN 

A-X 

ELSE 

B-X 

ENDIF 

V-tt 

FV-FW 

w-x 

FW-FX 

X-U 

FX-FU 

ELSE 

IF(U.LT.X)  THEN 

A-U 

ELSE 

B-U 
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END  IF 

IF(FU.LE.FW  .OR.  W.EQ.X)  THEN 

V-W 

FV-FW 

W-U 

FW-FU 

ELSE  IF(FU.LE.FV  .OR.  V.EQ.X  .OR.  V.EQ.W)  THEN 

V-U 

FV-FU 

END  IF 

END  IF 

11  CONTINUE 

PAUSE  ' Brent  exceed  maximum  iterations . ' 

3  XMIN-X 

BRENT-FX 

RETURN 

END 

C****  The  function  all  is  used  in  the  subroutine  BRENT  in  option  5.  *** 
FUNCTION  FUNC(P) 

COMMON  /XLIST/  PC,MC, F,MP, GAMMA, COV, VC, A,  TRAV, KINP, VEL, XBl, RHO, 
1  CVF,XLD,DVEL, IOPT,VCC,XMMC,D 

REAL  MC,MP,KINP 
MC*P 

VC-MC/RHO 
VC-VC/ (1.0-CVF) 

V1F-VC-C0V*MC 
VMF-VC+A*TRAV-COV*MC 
CONl-l+MC/ (2.0*MP) 

CON2-1+MC/(3.0*MP) 

CON3-GAMMA-1.0 

VBF-MC*F*CONl/ (GAMMA*C0N2*PC) +C0N3*V1F/GAMMA 
CON4- (VBF/VMF) **GAMMA 

KINP-PC/CON1*{GAMMA/CON3*VBF-V1F-VMF*CON4/CON3) 

VBFl-VBF+COV*MC-VC 

XBl-VBFl/A 

VEL-(2.0*KINP/MP) **.5 
FUNC— 1.0*VEL 
RETURN 
END 

Q*****************************************'*;******  ********************** 

C*********************  OUTPUT  SUBROUTINE  ****************************** 

Q*****1i*****1i****iiii**1t**1iiHi*iiii*li**ii*1i*1ilt**ii*1t***********1tii1i*****itit***** 

SUBROUTINE  XOUT 

COMMON  /XLIST/  PC,MC, F, MP, GAMMA, COV, VC, A, TBAV,K1NP, VEL, XBl, RHO, 
1  CVF,XLD,DVEL, IOPT,VCC,XMMC,D 

COMMON  /TIT/  TITLE, FILEIN,OUTFILE 
CHARACTER*20  TITLE*60,FILEIN,OUTFILE 
REAL  MP,MC,KINP 

Q******  ********  It*******  mi**  ***************  ******  *ii1iit1ti>1i1t***1i*1i***1iitit*1t* 

Q****************  First  Echo  of  Input  ********************************* 

Q* ********************************************************************* 

WRITE(18,799) 

799  FORMAT{//,T20,'  CONPRESS  PROGRAM  (Version  2)  RESULTS') 

WRITE (18, 800) TITLE 

800  FORMAT (//,'  Title:  ' , A60, //, 27X, '  INPUT  PARAMETERS',/) 

WRITE (18, 801)FILEIN,OUTFILE,IOPT 
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',A20,4X,' Output  File:  ',A20, 


801 

802 


803 

804 

805 


806 


807 

808 

809 


810 


811 


812 


813 


FORMAT ('  Input  File: 

1/,'  Run  Option:  ' ,12,/) 

WRITE (18, 802) 

FORMAT ('  Gun  Parameters' , 25X, 'Propellant  Parameters',/, 

1/ - ',25X,' - ',/) 

XVC-VCC/1.E3 

XF-F/1.E7 

WRITE (18, 803) XVC,XF 

FORMAT ('  Cham.  Vol.  (L) ', T23, ':' ,F10 .4, T35, ' Impetus  (J/g)', 
1T59,' :' ,F10.2) 

XD-D*10. 

WRI TE ( 1 8 , 8 0 4 ) XD , G AMMA 

FORMAT('  Tube  Diam.  (mm) ' , T23, ' : ' ,F10 . 4,T35, ' Gamma  (-)', 

1T59,' :',F10.5) 

WRITE (18, 805) TRAV, COV 

FORMAT('  Proj.  Trav.  (cm) ', T23, FIO . 4, T35, 'Covolume  (cc/g) ' 
1T59,' :',F10.5) 

XMP-MP/IOOO. 

WRITE (18, 806) XMP,RHO 

FORMAT('  Proj.  Mass  ()tg) ',  T23, ':' ,F10 . 5,  T35, '  Density  (g/cc) ' , 
1T59,' :',F10.5) 

XPC-PC/1.E7 

XMMC-XMMC/1000. 

WRITE (18, 807)XPC,XMMC 

FORMAT('  Max  Press.  (MPa) ', T23, FIO . 4, T35, '  Input  Prop.  ()cg) 
1T59, ' :' ,F10.5) 

XDVEL-DVEL/IOO. 

WRITE (18, 808) XDVEL, XLD 

FORMAT('  Target  Vel. (m/s) ',T23,' :', FIO. 4, T35, 'Target  LD  (g/cc) 
1T59,' :' ,F10.5,//,30X, 'RESULTS',/) 

WRITE (18, 809) 

FORMAT ('  Ballistic  Values' , 23X, 'Energy  Infomation' , /, 

1'  . ',23X,' . ',/) 

XV- VEL/ 100. 

XMC-MC/1000 . 

XPE-F/1 .E7*MC/ (GAMMA-1 . ) /I .E6 
WRITE (18, 810) XV, XPE 

FORMAT('  Velor  .ty  (m/s) ' , T23, ' : ' ,F10. 4, T35, ' Prop.  Energy  (MJ) ' 
1T59,' ,F10.5) 

XVC-VC/1000. 

XKINP-KINP*1.0E-07 

XPKE-XKINP/1.E6 

PER-XPKE/XPE*100. 

WRITE (18, 811) XB1,XPKE, PER 

FORMAT('  Burn  Out  (cm) ' ,T23, ' ,F10.4,T35,'Proj .  K£  (MJ) ' , 
1T59,' :' ,F10.5,'  (' , F5 . 1, ' %) ' ) 

XGKE-(1 ./3. )*MC/MP*XPKE 
PER-XGKE/XPE^IOO. 

WRITE (18, 812) XVC,XGKE, PER 

FOBMAT('  Cham.  Vol.  (L) ' ,T23, ' :',F10.5,T35, 'Gas.  KE  (MJ) ' , 
1T59,' :' ,F10.5,'  (' ,  F5 . 1, '  %) ' ) 

XGI-XPE-XPKE-XGKE 

PER-XG1/XPE*100. 

WRITE (18, 813) XMC,XGI, PER 

FORMAT{'  Prop.  Mass  (kg) ' , T23, ' : ' ,F10 . 5,T35, 'Gas  Internal  (MJ) 
1T59,' ,F10.5,'  (' , F5 . 1, ' %) ' ) 

XLLD-MC/VC 
WRITE (18, 814) XLLD 
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814 


FORMAT ('  Prop.  LD  (g/cc) ' ,T23, ' : ' ,F10 .5) 

ER- (VC+A*TRAV) /VC 
WRITE (18, 815) ER 

815  FORMAT(//,'  Expansion  Ratio  (-) ' ,T23, ' : ' ,F10.5) 

XXV- (VC+A*TRAV) /I .E6-COV*XMC/1000 . 

XGI-XGI*1.E6 

PMEAN-XGI* (GAMMA-1 . ) /XXV 
PBASE-PMEAN/ (1 . +1 . /3 . *XMC/XMP) 

PBASE-PBASE/1.E6 
WRITE (18, 816) PHASE 

816  FORMAT ('  Muzzle  Exit  Press  (MPa) : ' , FIO . 2) 

WRITE (18, 817) XMC/XMP 

817  FORMAT('  Mass  Charge/Mass  Proj  (c/m  ratio):  ',F10.5) 
PE-XKINP/XPC/TRAV/A 

WRITE (18, 818) PE 

818  FORMAT ('  Piezometric  Efficiency:  ',F10.3) 

RETURN 

END 
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Naval  Surface  Warfare  Center 

ATTN:  C«leG33. 

T.  Doi-an 
J.  Copley 

Code  C30,  Guns  and  Munitions  Division 
Code  G30L  D.  Wilson 
Code  G32,  Guns  Systems  Division 
Code  E23,  Technical  Library 
Dahlgren.  VA  22‘t48-5fl00 

1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  270PI.  Mr.  Ed  Oun 
101  Straus  Avenue 
Indian  Head.  MD  20610 

1  Commander 

Naval  Surface  Weapon  Cciuer 
Indian  Head  Division 
ATTN:  Cods  3120.  Mr.  Robert  Raa 
101  Straus  Avenue 
Indian  Head,  MD  20610 

1  Contmandcf 

Nava)  Surface  Weapon  Center 

Indian  Head  Division 

ATTN:  ckIc  21W^1.  Mr.  Ron  Simnms 

101  Suaus  Avenue 

In  lian  Head,  MD  20640 
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2  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  6210, 

Sharon  Boyles 
Norberto  Almeyda 
101  Straus  Avenue 
Indian  Head,  MD  20640 

1  Naval  Air  Warfare  Center 

ATTN;  Code  3891,  Mr.  Chan  Price 
China  Lake,  CA  93555 

1  Commander 

Naval  Air  Warfare  Center 

ATTN;  Code  3891,  Ms.  Alice  Atwood 

China  Lake,  CA  93555 

1  Commander 

Naval  Weapons  Center 
ATIT’  Code  388,  C.  F.  Price 
Info  Science  Div 
China  Lake,  CA  93555-6001 

1  OLAC  PL/TSTL 

ATTN:  D.  Shiplett 
Edwards  AFB,  CA  93523-5000 

10  Central  Intelligence  Agency 

Office  of  Central  Reference  Dissemination 
Branch 

Room  GE47  HQS 
Washington,  DC  20502 

1  Central  Intelligence  Agency 

ATTN;  Joseph  E.  Backofen 
HQ  Room  5F^2 
Washington,  DC  20505 

5  Director 

Sandia  National  Laboratories 
ATTN:  T.  Hitchcock 

R.  Woodfin 
D.  Benson 

S.  Kempka 
R.  Beasley 

Advanced  Projects  Div  14 
Organization  9123 
Albuquerque,  NM  87185 


2  Director 

Los  Alamos  National  Labaatory 
ATTN;  B,  Kaswhia 
H.  Davis 

Los  Alamos,  NM  875J5 

1  Directoi 

Lawrence  Livermore  National  Laboratory 
ATTN:  M.  S.  L-355,  A.  Buckingham 
P.O.  Box  808 
Livermore.  CA  94550 

2  Director 

Sandia  National  Laboratories 
Combustion  Research  FacUity 
ATTN:  R,  Armstrong 
S.  Vosen 
Division  8357 

Livennore,  CA  94551-0469 

1  University  of  Illinois 

Dept,  of  Mech./IndusL  Engr. 

ATTN:  Professor  Herman  Kricr,  144  MEB 
1206  N.  Green  St. 

Urbana,IL  61801 

1  The  Johns  Hopkins  Universily/CPIA 
ATTN:  T.  Christian 

10630  Little  Patuxent  Parkway,  Suite  202 
Columbia,  MD  21044-3200 

2  Pennsylvania  State  University 
Dept  of  Mechanical  Engr. 

ATTN:  Jeff  Brown 

312  Mechanical  Engineering  Bldg. 
University  Park.  PA  16802 

1  North  Carolina  State  University 
ATTN;  John  G.  Gilligan 

Box  7909 

1110  Burlington  Engineering  Labs 
Raleigh.  NC  27695-7909 

2  Institute  for  Advanced  Studies 
ATTN:  Dr.  H.  Fair 

Dr.  T.  Kichne 
4030-2  West  Baker  Lane 
AusUmTX  78759-5329 
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1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravenswood  Ave. 

Menlo  Park,  CA  94025 

1  SPARTA 

ATTN:  Dr.  Michael  Holland 
9455  Towne  Center  Dr. 

San  Diego.  CA  92121-1964 

5  FMC  Corporation 

ATTN:  Mr.  G.  Johnson 
Mr.  M.  Scale 
Dr.  A.  Giovanelti 
ivir.  J.  Dyvik 
Dr.  D.  Cook 
4800  East  River  Rd. 

Minneapolis,  MN  55421-1498 

2  Hercules  Inc. 

Radford  Army  Ammunition  Plant 
Manager  Manufacturing  Engineering 
Department 

ATTN:  D.  A.  Worrell 
Edvard  Sanford 
P.O.  Box  1 
Radford.  VA  24141 

1  Hercules  Inc. 

ATTN:  Dr.  Richard  Cartwright 
100  Howard  Blvd. 

Kcnvil.NJ  07847 

3  GT  Devices 

ATTN:  Dr.  S.  Goldstein 
Dr.  R.  J.  Greig 
Dr.  N.  Winsor 

5705A  General  Washington  Dr. 
Alexandria.  VA  22312 

3  General  Dynamics  Land  Systems 

ATTN:  Dr.  B.  VanDeu.scn 
Mr.  F.  Lunsford 
Dr.  M.  Weidner 
P.O.  Box  2074 
Wanen,  Ml  48090-2074 
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2  Alliant  Techsystems,  Inc. 

ATTN:  R.  E.  Tompkins 
J.  Kennedy 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

4  Olin  Ordnance 

ATTN:  V.  McDonald,  Library 
Hugh  McElroy 
Mr.  Thomas  Bourgeois 
Mr.  Dennis  Worthington 
P.O.  Box  222 
Sl  Marks.  FL  32355 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.  S.  Gough 
1048  South  St. 

Portsmouth,  NH  03801-5423 

1  Physics  International  Library 
ATTN:  H.  Wayne  Wampler 
P.O.  Box  5010 

San  Leandro.  CA  94577-0599 

2  Rockwell  International 
Rocketdyne  Division 
ATTN:  BA08. 

J.  E.  Flanagan 
J.  Gray 
6633  Canoga  Ave. 

Canoga  PaA.  CA  91304 

2  Princeton  Combustion  Research 
Laboratories,  Inc. 

ATTN:  M.  Summerfield 
N.  Messina 

Princeton  Corporate  Plaza 
1 1  Deerpark  Drive 
Bldg.  IV.  Suite  1 19 
Monmouth  Junction,  NJ  08852 

2  Science  Applications,  Inc. 

ATTN:  J.  Battch 

L.  TlionJhill 
1519  Johnson  Ferry  Rd. 

Suite  300 

Marietta,  GA  3O062.<>438 
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1  Eli  Freedman  &  Associates 

ATTN:  E.  Freedman 
2411  Diana  Rd. 

Baltimore,  MD  21209 

1  Rocketdyne 

ATTN:  Mr,  Otto  Heiney 
Mail  Stop  BA26 
6633  Canoga  Avenue 
Conoga  Park,  CA  91304 

1  Thiokoi  (Longhorn  Division) 

ATTN:  Dr.  David  Dillehay 
Mail  Stop  703-11 
P.O.  Box  1149 
Marshall,  TX  75671 

1  Thiokoi  (Elkton  Division) 

AT1T4;  Dr.  Rodney  Wilier 
55  Thiokoi  Road 
Elkton.  MD  21922 

1  Veriuiy  Technology,  Inc. 

ATTN:  Mr.  E.  Fisher 
4845  Millcrsport  Highway 
East  Amherst.  NY  14051-0305 

I  Veritay  Technology,  Inc. 

4845  Millcrsport  Hwy. 

P.O.  Box  305 

East  Amherst,  NY  14051-0305 

1  Baticllc 

ATTN;  TACTEC  Library.  J.  N.  Huggins 
505  King  Avc. 

Columbus.  OH  43201-2693 

2  California  Institute  of  Tcchnrriogy 
Jet  Propulsiort  Laboratory 
A1TN:  L.  D.  Strand.  MS  125-224 

D.  Elliot 

4800  Oak  Grove  Dr. 

Pasadena,  CA  91109 

I  General  Electric  Co. 

Defense  Sy.siwns  Divi.sion 
ATTN:  Dr,  J.  Mandzy 
Mail  Drop  43-220 
ICX)  Pla.stics  Avc. 

Pittsfield,  MA  01201 
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2  SAIC 

ATTN:  Mr.  N.  Sinha 
Dr.  S.  Dash 

501  Office  Center  Drive 
Fort  Washington.  PA  19034-3211 

1  State  University  of  New  York 

Department  of  Electrical  Engineering 
ATTN:  Dr.  W.  J.  Sargeant 
Bonner  Hall  -  Room  312 
Buffalo.  NY  14260 

1  Science  Application  International  Corporation 

ATIN:  Dr.  George  Chryssamellis 
8400  Normandele  Blvd. 

Suite  939 

Minneapolis,  MN  55437 

1  Sandia  National  Laboratories 

ATTN:  Mr.  Mark  Grubelich,  DIV  2515 
P.O.  Box  5800 
Albuquerque,  NM  87185 

1  IMl  Services  USA 

ATTN:  Mr.G,  Rashba 
2  Wisconsin  Circle 
Suite  420 

Cltevy  Chase.  MD  20815 


Aberdeen  Proving  Ground 

4  Cdr.  USACSTA 
ATTN:  S.  Walton 
0.  Rice 
D.  Lacey 
C.  Hcnid 

I  Dir.  USAHEL 
ATTN:  J.  Weisz 

6  Dir.  USARL 

ATTN:  AMSRL-WT-P.  Ingo  W.  May 

AMSRL-\fcT-T,  Walter  F.  Monisoo 
AMSRL-WT-PC. 

Ron  Anderson 
Dick  Beyer 
Joe  Hcimcrl 
Tony  Kotlar 
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54  Dir,  US.\RL 

AMSRL-WT-PA, 

William  Oberle  (25  copies) 

Avi  Birk 

Terry  Coffee 

Jim  DeSpirito 

Dawn  Gordner 

Arpad  Juhasz 

John  Kn^ton 

Charles  Leverilt 

Mike  McQuaid 

Kevin  Nekula 

Sharon  Richardson 

Irvin  Stobie 

Phuong  Tran 

Kevin  While 

Gloria  Wren 

AMSRL-WT-PB.  Ed  Schmidt 
AMSRL-WT-PD.  Bruce  Bums 
AMSRL-WT-PE. 

Andy  Brant 
I-ang-hiann  Chang 
Joe  Colburn 
Paul  Conroy 
A1  Horst 
George  Keller 
Doug  Kookcr 
Dave  Kmezynsky 
Rob  Lieb 
Tom  Minor 
Mike  Nusca 
Fred  Robbins 
Todd  Rosenberger 
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2  RARDE 
GS2  Division 
Building  R31 
ATTN:  Dr.  C.  Woodley 
Dr.  G.  Cook 
Fort  Halstead 

Sevcnoaks,  Kent  TN14  7BP 
ENGLAND 

1  Materials  Research  Laboratory 

Salisbury  Branch 

ATTN:  Anna  Wildegger  Gaissmaier 
Explosives  Ordnance  Division 
Salisbury  South  Australia  S108 

1  Laboratorio  Quimico  Central  de  Aimamcnto 

ATTN:  Capitan  Juan  F.  Hernandez  Tamayo 
Apartado  Oftcial  I.IOS 
28080-Madrid-SPAlN 

1  R&D  Department 

ATTN:  Dr.  Pierre  Aichambault 

S  Montee  des  Arsenaux 

Lc  Gardcur,  Quebec,  Canada  JSZ  2P4 

1  TZN  Fdrschungs-  und  Entwicklungszentium 

ATTN;  Mr.  T.  Wcisc,  Df.-lng. 

Abl.  SE,  Fachbcwich: 
Hochloislungspulsiechnik 
Ncucn.sothricthcr  StraBc  20 
D-SIO".  I  nu^lllb.  Ocmany 

1  Emsi'Mach-lnsuiul 

ATTN:  Df.  Gustav- Adolf  SchrOdcr 
Haupisurafic  18 

D-7858  Weil  am  Rheut,  Germany 

2  InsUiui  Ftanco-AHctnand 
ATTN:  Dr.  M.  Samirani 

Mr,  D.  Grutre 

F  68301  SAINT-LOUIS  C6dcx.  12. 
rue  dc  I'lndustric.  BP.  301.  France 
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